Functional renormalization for quantum phase transitions with non-relativistic bosons 
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Functional renormalization yields a simple unified description of bosons at zero temperature, in 
arbitrary space dimension d and for M complex fields. We concentrate on nonrelativistic bosons and 
an action with a linear time derivative. The ordered phase can be associated with a nonzero density 
of (quasi) particles n. The behavior of observables and correlation functions in the ordered phase 
depends crucially on the momentum kph, which is characteristic for a given experiment. For the 
dilute regime kph }t '^^'^ the quantum phase transition is simple, with the same "mean field" critical 
exponents for all d and M. On the other hand, the dense regime kph ^ ti}^'^ reveals a rather rich 
spectrum of features, depending on d and M. In this regime one observes for d < 3 a crossover to a 
relativistic action with second time derivatives. This admits order for d > 1, whereas d = 1 shows 
a behavior similar to the low temperature phase of the classical two-dimensional 0(2M)-models. 



I. INTRODUCTION 

The ground state of many quantum systems can undergo 
a second order phase transition if the density, concentra- 
tion, external fields or some effective coupling constants are 
varied. An example is the transition from para- to ferro- 
magnetism for bosonic atoms with spin. Bosonic quasi- 
particles also describe the quantum phase transition for 
strongly correlated electrons in case of anisotropic antifer- 
romagnetic order. We address here only systems where the 
low energy excitations are bosons - they can be relevant 
also for correlated electrons, if the fermionic excitations 
are gapped. 

An effective description accounts for a transition from an 
"ordered phase" with a nonzero continuously varying den- 
sity of bosonic excitations to a "disordered phase" where 
this density vanishes. At zero temperature the boson den- 
sity can often be characterized by a condensate, which is 
described by the nonvanishing expectation value of a (com- 
plex) field. In a broad sense such a condensate can be 
associated with order, while for the disordered phase the 
expectation value vanishes. Such transitions between a dis- 
ordered or "symmetric" phase and an ordered phase with 
spontaneous symmetry breaking are therefore described by 
Bose-Einstein condensation in a gas of interacting bosons. 
In this picture the parameter driving the phase transi- 
tion can be associated with an effective chemical poten- 
tial a. Universality of the critical behavior near second 
order phase transitions implies that many key features of 
quantum phase transitions are independent of the particu- 
lar "microscopic" physical systems. With ultracold bosonic 
atoms new ways of experimental investigation of such sys- 
tems open up. 

For dilute systems, many features of quantum phase 
transitions for bosons are well understood, and many de- 
tails of the critical behavior in various dimensions are 
known. The basic aspects are visible in a mean field theory 
and perturbation theory The use of several methods, 
including mapping to fermionic systems and bosonization 
in one dimension, together with strong universality argu- 
ments based on the renormalization group, allows for the 



computation of critical exponents and correlation functions 
. One finds a rather simple picture with mean field crit- 
ical exponents. This picture is valid, however, only as long 
as the system is sufficiently dilute. In this paper we ex- 
tend the discussion to dense systems. For one- and two- 
dimensional systems we will find qualitative changes. They 
are induced by the fluctuations of the Goldstone boson, 
which is characteristic for the spontaneous breaking of a 
continuous symmetry. In three dimensions, these effects 
are logarithmic - still sufficient for a cure of the infrared 
problems in many previous treatments This cure is 
similar to other renormalization group approaches [1, Q . 

In order to define the notion of dense and dilute, one 
should compare a typical physical length scale, I, with 
the average distance between particles, D ~ n~^l'^. For 
dense systems one has I :s> D, whereas dilute systems obey 
I <^ D. As a first attempt one could try to use the cor- 
relation length ^ as physical length scale. This works well 
in the disordered phase where ^ is finite away from the 
phase transition. In the ordered phase, however, the spon- 
taneous symmetry breaking of an abelian continuous sym- 
metry leads to superfluidity. For a nonzero condensate the 
system has always a gapless ( "massless" ) mode - the Gold- 
stone boson. The correlation length is infinite such that 
the system would appear "dense" for any nonvanishing n. 

In practice, a given experiment will always involve an 
effective momentum scale kph, for example the inverse of 
the wave length used to probe the system. Technically, 
the "physical momentum scale" kph may correspond to the 
momentum in some relevant Green's functions and act as 
an (additional) infrared cutoff for the fluctuations. The 
smallest possible value of kph is given by the inverse size of 
the experimental probe. Instead of we may therefore 
compare the physical momentum scale kph with the scale 
kp ~ n^/"^. For kp ^ kph the particle density only induces 
small corrections and the Bose gas is dilute. In contrast, 
the dense regime for kph kp corresponds to a situation 
where a characteristic inter-particle distance Z? ~ kp^ is 
small as compared to a typical experimental length scale 
I ^ k~^ . We use the concepts "dilute" and "dense" here in 
a rather general sense since we do not specify / a priori. 
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In the ordered phase the quantum fluctuations with low 
momenta, — , k kp, are dominated by the Gold- 
stone boson. We will call this the "Goldstone regime" . In 
three dimensions the Goldstone fluctuations play a quan- 
titative role, but do not change the qualitative behavior, 
except for the extreme infrared. (Due to a logarithmic 
running of dimensionless couplings the qualitative changes 
may only occur for exponentially small momentum scales.) 
Typically, fluctuations on length scales larger than the scat- 
tering length a give only small corrections. For kp^ 
the precise value of kpf^ becomes unimportant - the most 
important effective infrared cutoff is set by a~^. In this 
case we may consider a system with akp <C 1 as dilute, 
independent of kph- In contrast, for d = 1,2 the Goldstone 
fluctuations always play an important role for kph kp- 
In this case the value of kph matters and needs to be con- 
sidered as a separate physical scale. 

The simple critical behavior of the dilute regime always 
applies for the disordered phase since n — 0. In contrast, 
the understanding of the ordered phase is more subtle, in 
particular for lower dimensional systems, d — 1 or d — 2. 
The dilute regime kph 3> kp remains simple, with similar 
properties as for the disordered phase. As kph becomes 
smaller than kp we have to deal with a dense Bose gas 
where n sets a new scale. Following the scale dependence of 
"running" renormalized couplings one observes a crossover 
to the "Goldstone regime" , with new qualitative properties. 
We argue that for d = 1 and d = 2 the Goldstone regime is 
effectively described by a relativistic action with two time 
derivatives. It therefore shares common features with the 
classical 0(2Af)-models in dimension d+1. In particular, 
for d = 1 and M = 1 this implies the characteristic behav- 
ior of the low temperature phase in the Kosterliz-Thouless 
Q phase transition. 

We propose here a simple uniflcd picture for the prop- 
erties of the quantum phase transition which is valid on 
all scales. It is based on the functional renormalization 
group [7] for the average action Q , , ■ Within a sim- 
ple ^"'-model it describes the quantum phase transition for 
an arbitrary number of space dimensions d and an arbi- 
trary number of components M . Within the same model 
we can explore the flow in the disordered and the ordered 
phase. For low dimensions d = 1 or d = 2 we find several 
interesting crossover phenomena, indeed associated to the 
nontrivial physics of Goldstone bosons in low dimensions. 
This crossover persists for d = 3, but wide scale separations 
occur due to logarithmic running. The case d = 3 can be 
considered as the boundary dimension for the relevance of 
the Goldstone regime. 

All of the relevant physics is non-perturbative (with a 
few exceptions) and involves long range excitations. We 
do not limit our investigation to small interaction strength. 
We therefore rely heavily on the capability of modern ap- 
proaches to functional renormalization where the variation 
of an effective infrared cutoff enables the exploration of sys- 
tems with massless excitations (infinite correlation length) 
in a nonperturbative context for arbitrary d [^, [Tol. 
For the regime kph 3> kp the merits of our approach lie, for 
the time being, more in the simplicity of the unified pic- 



ture rather than in new quantitative results. In contrast, 
the fiow for the dense systems, kph <C kp, reveals features 
that have attracted less attention so far. 

Our approach is based on a functional integral formu- 
lation where the bosonic excitations are associated to a 
complex field x- It is formulated in a c? + I dimensional 
euclidean space with d space dimensions and an euclidean 
time T. (For nonzero temperature T euclidean time param- 
eterizes a torus with circumference T~^.) The transition 
from the Hamiltonian formulation with operators to the 
functional integral (or Lagrange formulation) with fields 
is sometimes subtle [Q, [ll], Two classes of systems 
can be distinguished, according to the presence of a linear 
T-derivative or not. In a rather general approach we may 
consider a microscopic or "classical" action 

S-Jx*{sd.-Vd^r-^-^)x + S.^, (I) 

X 

where / = J dr J d'^x. We will assume that Sint describes a 

local interaction, involving powers of x*X without deriva- 
tives. For arbitrary S and V the action ([T]) is invariant 
under euclidean time reversal r <-> — r , x ^ X*- 

The case 5 = is special, however. The system possesses 
now an enhanced rotation symmetry SO{d + 1), mixing 
space coordinates x and the time coordinate r. Indeed, 
a simple multiplicative rescaling of time or space coordi- 
nates brings the derivatives to the form [d"^ + A). The 
relativistic excitation spectrum can be directly seen by an- 
alytic continuation to "real time" , t = it. After suitable 
rescalings we may set V = 1, 2Mb = 1, such that eq. ^ 
reduces to the classical 0(2M)-model in dimension d+1, 
if X has M complex components. (For M > 1 suitable 
sums over components are implied in eq. ([T])). Functional 
renormalization has already provided a unified picture for 
the phase transition in classical 0(7V)-models for arbitrary 
d d, [13, mi , including the Kosterlitz-Thouless phase tran- 
sition for d -I- 1 = 2 and M = 1 [li]. Due to the en- 
hanced symmetry the vanishing of the coefficient linear in 
dr {S = 0) is stable under the renormalization fiow. 

In this paper we search for a similar unified picture for 
the "nonrelativistic bosons" with 5^ ^ 0. We will con- 
centrate on the simplest case V = where by a suit- 
able rescaling we may choose 5 = 1. In Minkowski 
space the microscopic action is now invariant under Galilei- 
transformations. We emphasize, however, that V — is 
not protected by a symmetry and second r-derivatives will 
be generated by the functional renormalization fiow. For 
T = this comes together with the higher order gradient 
terms requested by Galilei invariance. 

Even though we concentrate in this paper on the Galilei- 
invariant setting with V — 0, the discussion in sect. IVIIII 
will also cover the more general case of a microscopic ac- 
tion H]) with 5^0, V ^ 0. We first speciahze to one 
component and extend the discussion to M-components in 
sect. |Xl An overview over the different regimes for arbi- 
trary d, M and S/Vv, together with our main results, can 
be found in the conclusions. 
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Our paper is organized as follows. In order to specify our 
model and to fix notations we recall the functional integral 
formulation and functional renormalization in sect. HIl In 
sect, mil we derive the flow equations for the renormalized 
couplings for arbitrary d {M — 1). Sect. IIVI gives a brief 
description of running couplings in the disordered phase. 
We discuss the fixed points and the associated scaling be- 
havior relevant for the quantum phase transition in sect. 
IVl Sects. IIVI and fVl reproduce the known results for non- 
relativistic quantum phase transitions within the frame- 
work of functional renormalization. They may be skipped 
by the reader familiar with the subject. In sect. I VII we de- 
rive the flow equations for the ordered phase, first within a 
truncation where second r-derivatives are neglected. Sect. 
IVIII distinguishes the "linear regime" relevant for dilute 
systems from the Goldstone regime which is important for 
dense systems. For d < 2 our simplest truncation yields an 
attractive fixed point with nonzero order parameter and 
density. For d = I this fixed point will persist for extended 
truncations, while it turns out to be an artefact of the 
truncation for d > 1. 

For c? < 3 the large time behavior within the ordered 
phase is governed by a term quadratic in the r-derivatives 
~ V. We show in sect. IVIIII how this term is generated by 
the flow, even if it vanishes for the microscopic action. For 
dense systems and d < 3, the "relativistic dynamic term" 

V will always dominate over the term linear in dr- The 
running coupling S(k) vanishes with a power of k. In con- 
trast, for the boundary dimension d = 3 the vanishing of 
S(k) is only logarithmic. The corresponding change in the 
propagator induces qualitative changes for the renormaliza- 
tion flow and phase structure for one- and two-dimensional 
systems. They constitute the main result of the present 
paper. We find that the long distance behavior of corre- 
lation functions is similar as for the classical 0(2)-models 
in d 4- 1 dimensions. In particular, for d — 1 the ordered 
phase of the non-relativistic model behaves similar to the 
low temperature phase of the two-dimensional Kosterlitz- 
Thouless phase transition. We associate the Tomonaga- 
Luttinger liquid [l^ to this phase. In sect. |IX|we discuss 
the Goldstone regime in terms of non-linear cr-models. Our 
approach permits a unified view of the linear ip*-mode\s 
and non-linear cr-models. Sect. |X] is devoted to an exten- 
sion of our discussion to M complex bosonic fields. We 
present conclusions and outlook in sect. IXII 



II. FUNCTIONAL INTEGRAL AND 
FUNCTIONAL RENORMALIZATION 



We start with the partition for a nonrclativistic bosonic 
particle for M = 1 



Z 



r>xexp(-5[x]). 



(2) 



with action S[x] given by eq. (P), with S = 1,V = 0. The 
complex field x may be expressed by its Fourier modes 

Xix) = x{r,x) = J e^^^^xir,q) - J e^'^xiq), (3) 

with 



q •'q •'qo qo 



27r 



dqo, 



i2n)- 



d"- 



(4) 



For nonzero temperature T the Euclidean time r parame- 
terizes a circle with circumference VIt ~ T^^ and the Mat- 
subara frequencies go = iiinT , n e Z, are discrete, with 

Jqo ^ 

n 

In this paper we are interested in quantum phase transi- 
tions for T = 0. This phase transition occurs as the param- 
eter a is varied from positive to negative values. We regu- 
larize the theory by a momentum cutoff < and take 
A oo when appropriate. Furthermore, we assume the 
invariance of the classical action ([2|) under a global abelian 
symmetry of phase rotations x ~* e"''x, corresponding to 
a conserved total particle number 



N 



nix) = Vtd I n{q) 



(5) 



with the volume of d-dimensional space (fld+i = ^d^r)- 
Following the Nocther construction we can express n(q) by 
the two point correlation function 



n{q) = 



n, 



'+1 Jqo 



(6) 



We may associate a — a — Act with a chemical potential. 
Here our normalization of the additive shift Act [13| and of 
n{q) is such that iV = for ct < and TV 7^ for ct > 0. 
If we interprete n — N/Qd as the number density of some 
bosonic quasi particle, the quantum phase transition is a 
transition from a state with no particles to a state with 
nonzero particle density. 

By a suitable rescaling of units of x, t and x can 
replace 2Mb ^ 1 , ct ^ ct = 2MBO'/k^- Here we may 
use some arbitrary momentum unit k in order to make all 
quantities dimensionless, or we may retain dimensionful 
parameters by employing k = 1. (The parameters in Sint 
have to be rescaled accordingly, see Ref. [11].) In the fol- 
lowing we will work with a basis of real fields xi j X2 defined 
by x{x) = 7f (xi(a;) + iX2{x)) such that Xa{-q) = xliq)- 
The connected part of the two point function describes the 
propagator G 

{x*a{q)Xb{q')) - Qab{q,<^) + (x:(9)>(X6(g')>- (7) 

For a translation invariant setting, Q is diagonal in momen- 
tum space 



Qab{q, q') = Gab{q)S{q - q'), 



(8) 
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with 5{q - q') = (27r)'*+i(5(qo - 9o)^'^('? " 9")- Also, trans- 
lation invariance implies for a possible order parameter 
iXaiq)) = V2(f>o6{q)6ai with real 4>o- Here we have cho- 
sen the expectation value in the one-direction without loss 
of generality. 

We assume a repulsive two particle interaction (A > 0) 

S^nt = ^ljr{^)xi^)f ^^j^{Xa{x)Xa{x))\ (9) 

After the rescaling the mass dimensions are x ^ /i^^ , r ~ 
,_g ~ , 90 , X ~ > ~ /i^ , n - /i'', such 

that A ^ jjP'^'^. Already at this point one sees the crucial 
role of the dimension d. For d — ?> the coupling A has the 
dimension of a length. After a suitable renormalization it 
corresponds to the scattering length a ^ A. In the vac- 
uum (T = 0, n = 0, (7 = 0) the renormalized interaction 
strength A sets the only scale, besides the ultraviolet cut- 
off A. As a consequence, those correlation functions that 
are independent of A can only depend on dimensionless 
combinations, as \q and A^qo- For example, the two point 
function takes the form G — A{(f +^Sqo^-V q^)^^ with real 
functions A, S, and V/\^ depending on these dimensionless 
combinations. A nonzero density introduces an additional 
scale kp- The long distance physics will now depend on 
the dimensionless concentration c = akp ~ X'n}^^ . 

In one dimension (d = 1) the interaction strength scales 
A ~ /i. Now a length scale is set by A~^. For nonzero den- 
sity the macroscopic physics can depend on dimensionless 
combinations, as n/X. Indeed, for A — > oo the repulsion be- 
comes infinite such that a particle can never pass another 
particle. (This permits the mapping to a non-interacting 
Fermi gas 12;|.) The combination n/\ is a measure of how 
many interparticle distances a particle can travel before be- 
ing repulsed. It therefore defines an effective volume where 
it can move. 

The case d = 2 is special because A is dimensionless. 
In the vacuum no length scale except the cutoff is present. 
The running coupling vanishes logarithmically for large dis- 
tances (see sect. IIVI) . such that the long distance physics 
is described by a free theory. In the two point correlation 
G = A{(f + iSqo + Vq^)'^ the functions S, A can only de- 
pend on and A, besides cutoff effects involving q^/A^. 
No coupling V is allowed. A nonzero density sets again a 
further scale kp ^ n^/^. For momentum scales below kp 
the running of A gets modified and the macrophysics is 
no longer a free theory. In particular, the dimensionless 
combination Vn will now play an important role. 

We will conveniently work with the effective action r[0] 
which generates the IPI correlation functions. It ob- 
tains by introducing local linear sources j{x) for x(a;) and 
performing a Legendre transform of lnZ[j], with 4>{x) = 
{x{x))\j in the presence of sources 

m = -lnZ[j] (10) 

Jq 

The difference between S and F results from quantum 
fluctuations. We include these fluctuation effects stepwise 



by introducing first an infrared cutoff which suppresses the 
fluctuations with momenta q < k^. This is done by 
adding to the action ^ an infrared cutoff term Q 

AkS= I Rk{q)x*{q)x{q)- (11) 

Jq 

In turn, the effective action is now rep laced by the average 
action Tk which depends on fc Q, [ly|. With Rk{q) diverg- 
ing for fc — > 00 all fluctuations are suppressed in this limit 
and one finds Ffc^oc — S. On the other hand Rk{q) — 
for k ^ implies Ffc^o = T- The average action there- 
fore interpolates smoothly between the classical action for 
k ^ 00 and the effective action for k ^ 0. Its dependence 
on k obeys an exact flow equation Q 

dkTS] = iTr{afc7^(^f [^] + n)-^} (12) 

with TZ{q,q') = Rk{q)S{q — q'). The second functional 
(2) 

derivative F^ is given by the full inverse propagator in 
the presence of "background fields" 4>. For a homogeneous 
background field one has F^ {q,q') = P{q)5{q — q') with 
P a matrix in the space of fields (0i,02)- The trace in- 
volves a momentum integration and a trace over internal 
indices. Taking functional derivatives of eq. yields the 
flow of all IPI- vertices or associated Green's functions. Eq. 
p2p therefore describes infinitely many running couplings. 
For homogeneous background fields eq. takes the 

explicit form (with tr the internal trace) 

dkTk - %iir / dkRk{q){P{q) + Rk{q)y\ (13) 

The precise shape of the cutoff function Rk is, in principle, 
arbitrary. 

For dkRk decaying sufficiently fast for large (f the q- 
integration on the r.h.s. of the flow equation ([72]) or (|13p is 
ultraviolet finite. Instead of an explicite ultraviolet cutoff 
for the momentum integration, we can therefore define our 
model by specifying the form of F^, at some cutoff scale 

= A. The short distance physics is now given by the 
"initial value" Fa. For example, the bare coupling A in 
the action can be replaced by a coupling Aa, given by the 
fourth derivative of Fa. This definition has the advantage 
that momentum integrals can always be performed over an 
infinite range. The relation between the action S (with 
momentum cutoff) and the microscopic effective action Fa 
(without momentum cutoff) can be established by a one 
loop calculation. (In particular, this absorbs the shift Act 
in the chemical potential which is generated by the transi- 
tion from a Hamiltonian formalism to the functional inte- 
gral [13!] •) 

Our task will be to follow the flow of F^ from an ini- 
tial value given at fc = A towards fc = 0. From Vk=G = F 
the IPI-correlation functions of the quantum theory can be 
extracted by simple functional differentiation. Despite its 
conceptually simple one loop form, the exact flow equation 
(|12p remains a complicated functional differential equation. 
For approximate solutions we truncate the most general 
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form of Ffe. In the present investigation we will use very 
simple truncations, involving only a small number of cou- 
plings. The minimal set involves only three fc-dependent 
renormalized couplings S, rr? and A, according to the trun- 
cation 

Tfc = ^ ^4>*{Sdr - A + m2)0 + . (14) 

Nevertheless, many characteristic properties of the quan- 
tum phase transition in arbitrary dimension d will be ac- 
counted for by this truncation. This also holds away from 
the phase transition for the disordered phase and for the 
dilute regime of the ordered phase. For the dense regime 
of the ordered phase in c? = 1,2 we should add a term con- 
taining a second r derivative as in eq. ((1]), such that the 
minimal set consists of four running coupling m^. A, S and 
V . The coupling V is also needed for d = 3 if one attempts 
quantitative accuracy or a correct description of the long 
distance asymptotics. 

In this paper we mainly concentrate on non-relativistic 
bosons with a linear r-derivative in the action. The initial 
value Fa is then given by eq. p4)) . with S'a = 1 , to^ = —a 
and A = Aa- (We will only briefly comment on the more 
general case where a second r-derivative ~ Va is added to 
the microscopic action.) The microscopic average action 
Fa defines the model, which has only two parameters in 
our case, namely the rescaled chemical potential g and the 
microscopic interaction strength Aa . This should be distin- 
guished from the truncation of F^ for fc < A. In principle, 
all couplings allowed by the symmetries will be generated 
by the flow. This holds even though Fa has only two pa- 
rameters. Restricting F^ to a finite number of couplings 
defines the approximation scheme. 

In the formal setting the physical rt-point functions are 
only recovered for fc ^ 0. Nevertheless, the properties of 
Ffe for fc > also admit a physical interpretation. A typi- 
cal experimental situation has neither infinite volume nor 
observation devices working at infinite wavelength. This 
induces a characteristic experimental or "physical" momen- 
tum scale fcpTi, as mentioned in the introduction. Formally, 
this scale appears in the form of nonvanishing "external" 
momenta for the Green's functions which are relevant for 
a given observation. Often fcp/i acts as an effective infrared 
cutoff such that the evolution of these Green's functions 
(with finite momenta) stops once fc becomes smaller than 
fcp/i- On the other hand, for fc 3> fcph the external momenta 
are not relevant so that one may investigate the Green's 
functions or appropriate derivatives at zero momentum. 
In a simplified approach we may therefore associate the 
Greens-functions derived from F^^o at finite physical mo- 
mentum \q\ K, kph with the Greens function extracted from 
Tkj,^ at zero external momentum. In this picture we simply 
should stop the fiow of F^ at the physical scale kph rather 
than considering the limit fc ^ 0. The experimentally rele- 
vant Green's functions can then be extracted from F^^^ . Of 
course, such a procedure gives only a rough idea. In gen- 
eral, the Greens-functions will depend on several momenta. 
Even if only one momentum q is involved, the precise way 
how the flow is stopped by a physical infrared cutoff in- 



volves "threshold effects" [T3|- As a consequence, the pro- 
portionality coefficient between \q\ and kph will depend on 
the particular definition of the n-point function. 

The average action F^ has the same symmetries as the 
microscopic action, provided one chooses a cutoff Rk con- 
sistent with the symmetries. We sketch in appendix G the 
consequences of Galilei symmetry and local U{1) symme- 
try for the general form of F^ at T = - more details can 
be found in Ref. [l7|- 

III. FLOW EQUATIONS FOR POINTLIKE 
INTERACTIONS 

We first truncate the average action in the pointlike 
approximation and keep only the lowest time and space 
derivatives 

Ffc = f {Z^^*dr^-A4>*A^ + uiA^*^)} 

= j {S'0*9^0-0*A0 + u(0*0)}. (15) 

Here we have introduced S{k) — Z^[k)/ A{k) and the 
renormalized field 

(/) = A^/^^. (16) 

We use notations where quantities with a bar denote the 
couplings of the "unrcnormalizcd field" 0, whereas the 
renormalized couplings of the field cj) have no bar. At the 
scale k — A one has A{A) = Aa = 1 so that (f> and cj) co- 
incide. In general, the couplings and A are evaluated 
at a nonzero value of the renormalized field (pQ. As a con- 
sequence, they need not to be equal in the ordered phase, 
even in presence of Galilei symmetry for T = (cf. app. 
G). 

We choose the infrared cutoff function [l^ 

Rk = A{k^^q'')e{k'-q^). (17) 

This cutoff violates Galilei symmetry, but our truncation 
will neglect counterterms associated to anomalous Ward 
identities - they vanish for fc — > 0. The initial values of Fa 
will be taken as 

^0,A = ^A = 1 , UA = TOa0*0 -I- iAA(0*0)^; 

m\ = -CT. (18) 

Besides the rescaled chemical potential a our model de- 
pends on the strength of the repulsive interaction, Aa > 0. 
By a rescaling of the momentum unit k ^ k/a the parame- 
ters and fields scale as mj^ a'^niA , Aa a^^'^X^ , 
a^(j). Physical results for dimensionless quantities can 
therefore only depend on scaling invariant combinations 
as Aa(to^)~2~, q ^ /m\,qo/m,\,q ^ /A"^. We notice again 
the special role of d = 2 where Aa is dimensionless. 

The phase is determined by the properties of the effec- 
tive potential u for k —^ (or fc — > 1/L with L the macro- 
scopic size of the experimental probe) . In the ordered phase 
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the minimum of u occurs for (f>o{k) ^ and one observes 
spontaneous symmetry breaking (SSB) of the global C/(l)- 
symmetry. In contrast, the disordered or symmetric phase 
(SYM) has 00 = 0- For d = 1 we will encounter the bound- 
ary case where 0o(fc) ^ for arbitrarily small fc, while 
0o(fc = 0) = 0. (Typically (t>o{k) vanishes with some power 
of fc.) Since many properties of this phase are analogous 
to the SSB phase for d > 1 we will use the name "ordered 
phase" also for this case, even though long range order does 
not exist in a strict sense for the infinite volume limit. 

The flow of the average potential u follows by evaluat- 
ing eq. for space- and time-independent 0, i.e. eq. 
p3)) . We use the fact that the potential depends only on 
the invariant p — (jfcj) and Uk{p) = Vk{4') /^d+i- In our 
truncation the flow equation reads, using t ~ ln(fc/A), 



Here the propagator is a 2 x 2 matrix G ^ A ^G, 
+u' + 2pu" , -Sqo 



G-' - 



(19) 



(20) 



with q-' ^ q ^ ioT q ^ > and q'^ — k'^ for q < k^ . 
Primes denote derivatives with respect to p. Introducing 
the anomalous dimension 



Tj — —dt In A 



(21) 



we compute in app. A the flow equation for the average 
potential (at fixed (j) instead of fixed (j)) as 



dtu = ripu' + — -7/c 
ail 



1 - 



rf + 2 



k'^ + u' + pu" 



where 



Vfc^ + u'^k"^ + u' + 2pu" ' 



(22) 



(23) 



Eq. ([221) is a nonlinear differential equation for a function 
of two variables u[p, fc), if 'q{k) and S{k) are known. One 
may solve equations of this type numerically (ISy . 

We will choose here an even more drastic truncation and 
use a polynomial expansion around the minimum of u. In 
the symmetric regime the minimum of u is at p = and 
we approximate 



2 ^ A ? 

u = m p + 2^ 



(24) 



The corresponding flow equations for — m'(0) , A 
u"{0) read 



dtm^ 



rjm 



(25) 



dtX = 2r;A -f ^ ( 1 



V 



d+2j {k 



f^d+2 



Inspection of eqs. (jA.Sp . (|A.4p shows that the system ((25)) 
is closed and does not involve higher derivatives of the po- 
tential as u^^^ and u'^^K For the SSB regime, with minimum 
of u(p) at po 0, one expands 



(26) 



In this case the flow equations for po and A also involve u^^-* 
and u^^\ Neglecting these higher order couplings in our 
simplest truncation one flnds, from u'{po) — , u"(po) = 
A, the flow of the minimum 



dtPo = -jdtu'{po) 



(27) 



1 



y/k^ + 2Xpo + 2A/90 

The flow of the quartic coupling obeys now 

2vd ( V \ 9 k'^^^ 

' dS \ d+2 ^WT2Xio 



(28) 
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1 3 ^\ f ^ 

In the symmetric regime we flnd a (partial) flxed point 
for m? = 0, while in the SSB one has a flxed point for 
Po = 0. These points coincide, with a quartic potential 
u — \Xp^ . In turn, the flow for the quartic coupling 



(29) 



has a flxed point for A = 0, corresponding to a free the- 
ory. In order to understand the flow pattern we will need, 
however, the flow of appropriately rescaled dimensionless 
quantities and the behavior of S and 77. 

For a computation of rj and 5* we need the flow of the 
inverse propagator 



ah 



S<j,l{q')SMq") 



Pab{q')S{q' - q" 



(30) 

The flow of Pab obtains by the second functional derivative 
of the exact flow equation (fT^ 

dtPabiq) = / dtRk{q'){G\d{q') (31) 



with 



iladelbfcGefiq' + q) + Ibdelaf cG ef {q' - q)} 



G{q') = {P{q') + Ru{q')) \ 



(32) 

We have omitted here a term ^ F^''^ which docs not con- 
tribute to 77 or 5* in our truncation of momentum indepen- 
dent vertices. The cubic couplings (j)^ are specified by 



^rfi(y,y') 
Hl{q') 

5<t>b{q") 



7acd(l)S{p - p" + q), 

ibcdW - p" - q") (33) 
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and read in our truncation 



Jacd = 7acd = V2A'^^u"{5alScd + SdSbd + SdlSac) 

+2pu'^''^Sai5ci5di}. (34) 

The anomalous dimension t] and the flow of S are defined 
by 



and 



V = 



-dtlnS 



1 d 



dtP22{q)\q=o 



1 d 
SA dqo 



dtP2l{q)\ 



(35) 



(36) 



Many quafitative features for arbitrary dimension d can 
already be seen in the extremely simple truncation of this 
section. Nevertheless, an important missing ingredient for 
the dense regime is the second order r-derivative ~ V4)*dl;.(t> 
discussed in sect. IVIIII and appendix C. We will see that 
this plays a central role for the infrared behavior, and ap- 
propriate corrections should be included in the flow equa- 
tion ([22|) . Beyond this, the extension of the truncation 
is more a matter of quantitative improvement. The most 
general pointlike interactions are accounted for by eq. (|15p . 
For example, including in m a term ^ describes point- 
like six-point vertices, as discussed in app. E. The leading 
order in a systematic derivative expansion needs, in addi- 
tion to V , a term — iFpAp— jYtpd'^p with p = (l>*(f>. This 
contains momentum dependent interactions. The next to 
leading order in the derivative expansion has A, S, V and 
Y,Yt depending on p. All these approximations have been 
successfully implemented for "relativistic" models with sec- 
ond order T-derivatives and have led to a precise picture 
for 0(iV)-models in arbitrary d [ISl, [ll. 



IV. DISORDERED PHASE 

In the next four sections we discuss the simplest trunca- 
tion. In the symmetric regime ((/)o — 0) the cubic couplings 
~ (j)"f vanish. From eq. (|3ip we find in our truncation of 
momentum independent vertices 



?7 = , T]s =0. 



In terms of the dimensionless mass term and quartic cou- 
pling 



= mVfc2 , A = 



Afc 



d-2 



(38) 



we obtain 



am^ = 



dtw = —2w, 



dt\ = {d-2)X + ^{l + wy^X^. (39) 

Since m? , A and S do not depend on fc, the quantum prop- 
agator (F^^^)"^ is given by the classical propagator (for real 
frequencies oj — —iqa) 



G = {-uj- 



(40) 



The non-renormalization property of G for T = 0, > 
is believed to be exact since the situation describes the 
vacuum with zero particle number 12^. This is also the 
reason for the closed form of eq. (|39p which docs not involve 
higher order n-point functions. 

As long as fc^ ^ (or w <C 1) the quartic coupling 
A runs while for k'^ <C m?{w 3> 1) the running effectively 
stops. For d <2 and w — Q the combination A is attracted 
towards an infrared fixed point at 



A* 



(2 - d)d 
4wrf 



(41) 



In the vicinity of this fixed point A decreases with k 

Xr^kk^-'^ (42) 

and the repulsive interaction tends to be shielded by the 
fluctuation effects. For d > 2 there is no fixed point for 
A 7^ 0. Again A decreases for k ^ 0. Now the running of A 
stops in the infrared even for w = 0. 

The explicite solution of the flow equation for A in the 
range w 1 reads for d 2 



X{k) = Aa 



4fdAA 
d{d~2)S 



(A 



d-2 



■d-2 



(43) 



We note the different behavior for c? > 2 and d < 2. For 
d > 2 the fluctuation effects on A are dominated by the 
short distance physics, i.e. momenta of the order A (ultra- 
violet domination). One expects the precise value of the 
effective quartic coupling to depend sensitively on the mi- 
croscopic details. In contrast, for d < 2 the long-distance 
physics dominates (infrared domination) . For systems with 
a characteristic physical infrared cutoff kph the value of the 
effective coupling is given by X{kph)- For d <2 the correc- 
tions are dominated by the fluctuation effects with infrared 
momenta « fc^^ . If the microscopic coupling Aa is large 
enough, Aa > Xc{kph), 



^37) the value of Aa becomes unimportant 



X{kph) ~ Xc{kph). 



(44) 



(45) 



The system has lost memory of the microscopic details ex- 
cept for the value of m\. For w = the value Xc{kph) is ac- 
tually an upper bound for the allowed values of X{kph). For 
kph ~^ the model becomes non-interacting, X{kpii) 0. 
This "triviality property" is analogous to the relativistic 
model, as relevant for the upper bound on the Higgs mass 
in the standard model of particle physics. For > one 
effectively replaces k^f^ cnP' with c a proportionality 
constant of order one. 

The boundary between the qualitatively different role 
of fluctuations occurs at the "upper critical dimension" 
dc — 2. For d > dc the critical behavior is well approxi- 
mated by mean field theory, with mean field theory critical 
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exponents. For d < dc the fixed point behavior (|4T|) influ- 
ences the critical physics as far as the interaction strength 
is concerned. At the critical dimension d ~ 2 the running 
of A for w ^ 1 becomes logarithmic 



A(fc) = Aa 



Aa A' 
AttS k 



(46) 



V. 



SCALING SOLUTIONS AND QUANTUM 
PHASE TRANSITION 



It is instructive to discuss the critical behavior in terms of 
the scaling solutions. Possible scaling solutions correspond 
to the fixed points for w and A, i.e. to values where both 
dfW and dtX vanish. For all d one has the trivial fixed point 



(A) : = , A* = 0. 



(47) 



Small deviations from this fixed point, with w > 0, grow for 
fc — > 0. The fixed point (A) is unstable in the w-direction, 
thus w (or m^) is a relevant parameter. For d > 2 fixed 
point (A) is infrared stable in the A-direction, i.e. A is an 
irrelevant coupling. However, for d < 2 A the coupling 
becomes a relevant parameter, too. Fixed point (A) has 
two IR-unstable directions for d < 2. Indeed, the flow of A 
is attracted towards a second fixed point 



(B) 



= , A = A, 



(48) 



with A* given by eq. (|¥T|l . Fixed point (B) has only one 
relevant parameter w whereas A becomes irrelevant. The 
critical behavior is dominated by fixed point (B), except for 
very small Aa where one observes a "crossover" of the flow 
from the vicinity of (A) towards (B). Both fixed points 
(A) and (B) are located exactly on the phase transition 
rn^ = 0. 

The value of A does not affect the flow of w or the anoma- 
lous dimension t] or ris- We therefore find for the symmet- 
ric phase a mean field critical behavior for m\ —> 0. This 
equally applies for both fixed points (A) and {B) which are 
distinguished only by the value of A. There is no running 
of and the anomalous dimension rj as well as r]s vanish. 
The correlation length ^ = m^^(k — > 0) simply obeys 



— ^ = 1^1- (49) 



m{k —> 0) ?7iA 
and the correlation time (for m\ > 0) is given by 



(50) 



The time averaged correlation function for mj^ — decays 
according to the canonical dimension {d > 2) 

{^*ir)m)-M'~', (51) 

as given by the d-dimensional Fourier-transform of eq. (j40p 
for m\ = 0, a; = 0. 



The corresponding critical exponents are the mean field 
exponents [2l[, 0] 



77 = 



(52) 



In the present case, the critical exponents follow from naive 
scaling arguments. More generally, the critical exponent 
ry corresponds to the anomalous dimension for the scal- 
ing solution. Indeed, if we evaluate the propagator for 
(f ^ > 0, the external momentum acts like an infrared cut- 
off {\q\ ~ kph), such that A ^ k-'> ^ (g 2)-'7/2. At the 
phase transition the static propagator (go = 0) behaves as 
G = G/A ^ {q '^)~^~^^^'^ , which is precisely the definition 
of the critical exponent ry. 

The value of rjs for the scaling solution determines the 
dynamical critical exponent z. 



(53) 



The dynamical critical exponent z relates the <f- 
dependence and the go-dependence of the renormalized in- 
verse propagator away from the phase transition 



G~\qo = 0,q)=q^ + r\ 
G-\qo,q^O)^lS{qo)qo+r^ 



(54) 



icq^' 



If the zeros of G~^(go) occur for a value of qo with posi- 
tive real part, i?e(qo) — '''^^^ the correlation function for 
real time t decays exponentially with a typically dissipation 
time Tc, implying for <f ^ <C ■^^^ 



{ip{t,q)ip*{0,q)) -exp(-i/rc). 



(55) 



Assuming that for the zero of G~^ one has Re{qo) ^ Im{qo) 
one can relate the dissipation time Tc to the correlation 
length ^ 



(re) 



2/z 



(56) 



A nonzero external qo will replace the infrared cutoff in 
the propagator, fc^ — > S{qo)qQ, such that 



Siqo)qo - k-'^'qo ^ [5(go)go]"''"/'go. 

The scaling 



{S{qo)qo 



qo , S{qa) go 



(57) 



(58) 



yields the relation ([55]) between z and 775. A simpler ar- 
gument compares the scaling of a characteristic go with 
, go ~ fc^, where go is determined such that the go- 
dependent part in G^^ has the same size as the IR cutoff 



S{k)qo ^ e 



-"-go 



90 



(59) 



This yields, of course, the same relation (|53p. 

Obviously, the scaling arguments leading to eq. ([55]) 
depend strongly on the absence of any other relevant cut- 
off. They will not be valid for dense systems where \/2Apo 
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constitutes an important infrared cutoff for the radial fluc- 
tuations. For the dense systens we find for all d a small 
momentum behavior G"^ ^ {([^ + Qo/v'^), such that the 
same type of scaling arguments yields z = 1. 

For d > 2 mean field theory is expected to be a valid ap- 
proximation. For d < 2, however, the strong dependence of 
A on fc will result in a momentum dependence of the effec- 
tive vertex, with k replaced by \fW^ ^ ^^^d p a characteristic 
external momentum of the vertex. The approximation of 
a pointlike interaction becomes inaccurate and one may 
question the validity of a mean field tretament. Neverthe- 
less, relations dtva^ — 0,7] — r]s — continue to hold (cf. 
eq. (001)), implying the mean field critical exponents ([5^ 
for all d. Also the equation for a momentum dependent 
quartic coupling will remain closed. Only the value of A* 
and the precise evolution of the quartic coupling A will be 
modified for extended truncations. For d — 1 and Aa — > oo 
our model corresponds to "hard core bosons" . In d = 1 
this is equivalent to a model of free spinless fermions and 
the universality class for fixed point (B) is therefore known 
, confirming that eq. ((52)) is exact for d = 1 . 

The quantum phase transition at = is the only 
phase transition that we discuss explicitely in the present 
paper. Its scaling properties in the symmetric phase are 
quite simple. Our functional renormalization group equa- 
tions account well for these scaling properties, establishing 
them as a reasonable starting point for T > in a straight- 
forward generalization where the go-integration in the ap- 
pendices A and B is replaced by a Matubara sum. The 
simple features of the quantum phase transition discussed 
in this paper, namely the location exactly at n = 0, the 
non-renormalization of rri^ and the vanishing rj and rjs, are 
all particular for the non-relativistic models with V = 0. 
(We will see in I Villi that = is stable with respect to 
the fiow in the symmetric phase.) Starting with V ^ will 
change all these properties - for example, the relativistic 
model with V , S = shares none of them. We also do 
not investigate here the possibility that a very strong repul- 
sion annihilates the order even ior V — , T — , n 0. 
(This would lead to a different type of quantum phase tran- 
sition at some critical value of Aa, perhaps characterized by 
the critical behavior of the relativistic 0{2M)- or U{M)- 
models.) 

We close this section by a remark that a line of fixed 
points exist for all values m\ > 0. The associated scaling 
solutions reflect, however, a different scaling behavior 
Indeed, for fc^ <C we may use the variables and 



A: 



A. 



From eq. (I25|) and for f] = rjs = we extract 



dtX= id + 2)X+-j^X 



(60) 



(61) 



and observe that an infrared stable flxed point A = exists 
for all d > 0. As fc^ crosses the "threshold" rri^ the flow of 
A shows a crossover from fixed point (B) (or (A)) to the 
fixed point of eq. (pT|) . This is, of course, a fancy way of 
stating that A stops running. 



VI. ORDERED PHASE 

We next turn to the ordered phase. This will be charac- 
terized by a richer spectrum of physical phenomena, since 
even for T = the particle density is nonvanishing. We will 
see that for c? < 2 the long distance physics is always char- 
acterized by an effective theory with strong interactions. 
The quantum phase transition to the disordered phase re- 
mains simple for d = 3 and small coupling since fluctuation 
effects play a minor role. Such a simple description also ap- 
plies for d < 2 in the dilute regime, as long as the momenta 
and energies of the process considered are larger than a 
characteristic momentum kp ot a characteristic energy ep 
related to the density. (In our normalization ep = kp.) 
For smaller momenta and energies, however, the density n 
sets a new scale 



kp — 



dn \ 
8^d) 



l/d 



(62) 



(The normalization of kp ~ ni/'' is somewhat arbitrary and 
we have chosen it here in analogy with a Fermi gas of par- 
ticles with spin 1/2.) One expects that this scale strongly 
influences the long distance behavior. For d = 1, 2 one finds 
qualitatively new phenomena whenever the "physical mo- 
mentum" kpfi is smaller than kp (dense regime) . For d = 3 
and small couplings these effects may matter only on ex- 
ponentially small scales, since the running of dimensionless 
couplings is logarithmic. 

The new physics for non zero density is directly related 
to the possibility of a condensate. The flow equations 
will be influenced by the "renormalized order parameter" 
Poik) > 0, which denotes the value of p for which the av- 
erage potential Uk is minimal. For T = the symmetries 
relate the asymptotic value po — po{k = 0) to the density, 
cf. app. G, 



n = Po. 



(63) 



(For dimensionless fields this reads n = pok'^ .) In terms of 
the original fields the asymptotic order parameter po — 
Po{k = 0) denotes the condensate density ric 



n-c = Po, 

such that the condensate fraction ilr reads 



1), = ^ = ^ 
n Po 



A- 



(64) 



(65) 



We will encounter the notion of a "local condensate" po{k), 
even if no long range order exists, i.e. if (f>o{k —^ 0) = 0. 
In this perspective k can be associated with the inverse 
size of a domain and (j)o{k) measures the expectation value 
of the order parameter in such a domain. We refer to the 
"SSB- regime" of the flow whenever po(fc) 0- The ordered 
phase or the phase with spontaneous symmetry breaking 
(synonymous) are characterized by a nonzero po at the end 
of the running, i.e. for fc = or L~^. Technically, the 
running of the couplings in the SSB-regime is more involved 
due to the presence of the cubic couplings for (f>o ^ 0. 
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In the SSB-regime the two modes (pi and 4>2 show a dif- 
ferent behavior. With an expectation value (f>o in the 1- 
direction, 0i denotes the radial mode (longitudinal mode) 
which is typically "massive" or "gapped" , with "mass 
term" 2Xpo. In contrast, the "Goldstone mode" (transverse 
mode) 4)2 is massless. In the SSB-regime, the relative size 
of the contributions from the Goldstone and radial modes 
is governed by the dimensionless ratio 



w 



(66) 



In eq. ^7^, we note that for w =_2 the flow of the 
unrenormalized parameters po = p/A, X — A^X vanishes. 
For w < 2 (and t] < d+2) one finds for fc — ^ an increasing 
Po and increasing A, whereas for w > 2 both quantities 
decrease. 

Within the SSB-regime we distinguish between two lim- 
iting cases. The "linear regime" refers to w <C 1 where both 
the radial and Goldstone mode are equally important. In 
contrast, in the "Goldstone regime" for w ^ 1 the radial 
mode plays a subleading role and the dominant physics is 
related to the behavior of the Goldstone modes. 

From eqs. ([27]), (|25|) we extract the flow of w 



the dynamical critical exponent z in eq. (j53p holds only as 
long as the first order r-derivative dominates. For V ^ 
and ?7s < — 1 the relativistic dynamic term will domi- 
nate, yielding simply z — 1. Furthermore, in the ordered 
phase Goldstone bosons dominate the correlation function 
at large distances in space or time. We will see in sect. 
X that the decay of the correlation function is powerlike, 

~ [9*^ + {qo/v)"^]^ ^^^y leading effectively to z = 1, 
independently of all other details. We observe an apparent 
clash with the critical exponent z = 2 in eq. ((52)) . 

We may insert our results for 77 and rjs into eqs. (|67p . 



dtw 



dtX 




XU-2 + ^A(2 - wf{l + w)-''/^ 
a 

2vjX'^w 4 - 6w - I 



d^{d + 2) (l+wy 



(72) 



a J o , , ^^'^ ?7 \ Aw 

Otw = w < —2 + 77 H 1 , 

1 ' 2d \ d + 2 J y^TT^ 



1 - 



1 + w 



1 



1 + w 



(67) 



In eq. (j67p we encounter again the dimensionless combina- 
tion A = Xk'^^^/S. Its evolution obeys (77 = —dt In A, 775 = 
-dtlnS) 



dtX = {d~2 + 27^ + r]s)X 



1 



1 - 



2«d A ?/_ 

d \ d + 2) yrTw ^ + w 

l_2w-^~^^ 
l + w 



(68) 



The anomalous dimension 77 is computed in appendix B 
and we find in our truncation 



r^='^J!£\w{l + w)-^'\ 
d 



(69) 



It vanishes both for w ^ and w 00. For 775 we find 
(app. B) 



^ ^A77,^-i^^. (70) 



"^"^ 2dy d+2)''"' (i + w)5/2 

We note that the leading term for large w 



Vs 



2d V d+2 



(71) 



cancels in the flow of A the term ^ A^w^/^. 

We note that t;^ can take large negative values. In this 
context we emphasize that the relation between 775 and 



These two coupled nonlinear differential equations for the 
two couplings w and A already yield several characteristic 
features of the ordered phase for arbitrary d. However, the 
understanding of the dense regime requires an extension 
of the truncation by the "relativistic dynamic term" ~ V. 
This will be necessary in order to get the correct behavior 
for -u; ^ 00 (sect. IX). 



VII. LINEAR AND GOLDSTONE REGIMES 

It is instructive to consider first the linear regime in the 
limit iti — > 0, where 



dtW 
dtX 



-2w 



1 + ^Xw 
d 



^^^^ -(At.)^ 



d-2 



d^{d + 2) 



d2(d + 2) 



(73) 



For fc — > one finds that w increases. The dimensionless 
interaction strength A decreases for d > 2, while it increases 
for d < 2 and small A. For w = we recover the two 
fixed points (A) and (B) already found in the symmetric 
phase. As before, (A) is IR-stable in the A-direction for 
d > 2 and unstable for d < 2. Fixed point [B] at w^, — 
, A* = (2 — d)d/{Avd) exists for d < 2 and is IR-stable 
in the A-direction. For both fixed points is a relevant 
parameter. For fixed point (B) we find for k ^ that 
2Xpo — wk"^ — M^ approaches a constant, as well as S and 
A. Since A — > SX^k^~'^, we find that po(fc) increases for 
d < 2 and w <C 1 according to 



(74) 
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This behavior stops once w reaches a value of the order 
one. 

The quantum phase transition occurs for w — 0. The 
critical behavior is characterized by fixed point (A) for d > 
2 and {B) for d < 2. Since the fixed points are the same 
for the ordered and disordered phases we also obtain the 
same scaling behavior. At this point everything may look 
very simple. 

A closer look at the ordered phase reveals, however, that 
fixed point {A) or (B) cannot describe all aspects of the 
quantum phase transition. It is not clear how the expo- 
nents v, rj and rjs should be defined in the ordered phase. 
The correlation length ^ could be defined in the radial di- 
rection, = (2Apo)~^^^ where Apo should be evaluated for 
k = . In the Goldstone direction, however, the correla- 
tion length is infinite. We emphasize that the correlation 
function for the complex field is dominated by the prop- 
agator for Goldstone bosons G22 which does not exhibit a 
finite correlation length but rather shows a powerlike decay 
for large |x|, 



(0*(a;)<^(O))-^i 



(0i(a;)0i(O))-02^-(02(x)02(O)) 



(75) 



Since Gn decays faster than G22 only the latter matters 
for \x\ — > 00. 

Similarly, Goldstone bosons dominate the occupation 
number (l6l) for the small momentum modes 



n{q) 



{G22{qo,q) + Gii{qo,q)) 



1 

' 2' 
(76) 

We will see that the shape of G22 for small q ^ and go 
can become nontrivial and is no longer governed by the 
"quantum critical fixed point" (A) or (B). 

We may also study the critical behavior of the (bare) 
order parameter pQ{k 0) ^ a^^^. The fiow for fc ^ 
will necessarily involve the flow in the region of large w 
and one may wonder if this can be described by the fixed 
points (A) or (B) anymore. We will even find that for d — 1 
the order vanishes in a strict sense, po{k ^ 0) ^ 0. The 
definition of (3 seems not to be meaningful anymore. These 
particularities of the correlation length and the bare order 
parameter cannot be explained by simple extrapolations 
from the fixed points (A), (B) for which Goldstone bosons 
play no particular role. 

What actually happens is a crossover phenomenon be- 
tween the scaling associated to fixed points {A) or (B) for 
the quantum phase transition and the Goldstone regime 
where the gapless Goldstone mode dominates. This 
crossover depends on the scale kph of characteristic mo- 
menta of an experiment. The crucial quantity is the ratio 



w{kph) 



2X{kph)po{kph) 



(77) 



Only for 'w{kph) ^ 1 are the scaling laws of the quantum 
phase transition given by the fixed points [A) or {B). In 



the opposite limit one has to explore the Goldstone regime 
w ^ \. We discuss in appendix F that -y/2Apo plays the role 
of the momentum scale kp ^ n^^''' associated to the den- 
sity. The linear regime therefore applies for dilute systems, 
while the Goldstone regime is relevant for dense systems. 
We notice that for any small nonzero order parameter po 
(corresponding to a situation near the phase transition) 
there is always a range of very small momenta kph such 
that w(kph) 3> 1. The extreme long range behavior is al- 
ways dominated by the physics of Goldstone bosons for 
which the fixed points {A), [B) are not revelant. This also 
matters in practice since the macroscopic size of an exper- 
iment corresponds to very small kph- 

We show the different regimes which are relevant for the 
quantum phase transition in Fig. 1. The "linear" regime 
is governed by the fixed points [A) or {B). In the disor- 
dered phase this finds a simple extension to the "massive 
regime" where the fiow simply stops due to the presence of 
an infrared cutoff ~ m, cf. eq. (j6ip . In the ordered phase, 
however, the "Goldstone regime" is qualitatively different. 
The long distance physics shows new features which cannot 
be explained by the fixed points {A), (B). We note that 
for a fixed nonzero kph the immediate vicinity of the phase 
transition for cr ^ is always governed by the quantum 
critical point, while for fixed a > one always enters the 
Goldstone regime as kph 0. 
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FIG. 1: Schematic view of different regimes for quantum phase 
transitions. Dense systems are described by the Goldstone 
regime. 

We will next explore the Goldstone regime. This will 
bring us to the most important results of this paper. For 
the qualitative behavior of the ffow away from the exact 
location of the phase transition (p ^ 0) we observe an 
important difference between the disordered and ordered 
phase. In the disordered phase the running of the couplings 
rr? , A simply stops once fc^ 1 and the consequences of 
scaling are immediate. In the ordered phase, however, we 
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encounter the massless Goldstone fluctuations at all scales, 
including fc^^ ^ 1. Correspondingly, the flow equations in 
the regime w 3> 1 will be nontrivial and we should explore 
their consequences. 

We first work with our simplest truncation and extend 
it subsequently in the following sections. Within the trun- 
cation (jlSp. (pS)) we will find a new fixed point of eq. ([7^ 
for d < 2 and a nontrivial scaling behavior for 2 < d < 3. 
As more couplings are included in extended truncations we 
find that the fixed point persists for d = 1, while it turns 
out to be an artefact of the truncation for d > 1. The pre- 
cise properties of the fixed point are quite sensitive to the 
truncation, and the "lowest order results" of the simplest 
truncation have to be interpreted with care. Indeed, for 
< c? < 2 the flow equations in the ordered phase ([72]) 
exhibit an additional fixed point for w^, 

(C) : w,^Q,~X,^ 0. (78) 

The characteristic fixed point values obtained in this trun- 
cation by a numerical solution of eq. (j72p are shown in 
table I. In appendix E we have performed a similar compu- 
tation including a local six point vertex in the truncation. 
Comparison of tables I and II reveals a very strong trunca- 
tion dependence for d near two, while the results are more 
robust for d = 1 . 



d 






V 


Vs 


<* 


d + r]s + 1J 


1 


2.257 


28.3 


1.73 


-2.65 


-8 


0.08 


1.9 


7.33 


23.8 


0.335 


-2.012 


-80 


0.23 


1.99 


127.3 


5.85 


0.021 


-1.99 


-500 


0.021 



table I: Fixed point values for (C) 



Fixed point (C) is infrared attractive in all directions. 
Within the restricted space of couplings considered in our 
simple truncation this is an example of "self-organized crit- 
icality". For < d < 2 the flow for any initial value 
m\ < , Aa > will finally end in fixed point (C). In ta- 
ble 1 we also indicate a characteristic value oi t = ln(fc/A) 
for which the fixed point is reached (for initial w and A 
of the order one). As d approaches two the fixed point 
behavior sets in only at extremely large distances. 

Since for fixed point (C) w and A take constant values 
one finds in the simplest truncation 

Po{k) = , A(fc) = KS{k)k'-''. (79) 

With 
we infer 

For rjs < — d the rcnormalized order parameter po{k) in- 
creases with k, while for d + rjs + rj > the bare order 



parameter po = A~^pQ vanishes for fc — > 0. From the val- 
ues of d -|- 77s -|- 77 in tabic I one would infer that no long 
range order is present for d < 2. (For the one-dimensional 
boson gas we find in this simple truncation that po van- 
ishes ~ k'^-'^^.) Then there is no meaningful definition of 
the critical exponent /3 for d = 1. Also 

^j,{k)^{2X{k)po{k))-'^' = -L^ (82) 

always diverges for fc 0, due to the existence of the 
fixed point for d < 2. For such a behavior there would be 
no meaningful definition of a correlation length even for 
the radial mode, due to the strong impact of Goldstone 
fluctuations. 

Within our simplest truncation one would conclude that 
for d = 1 a quantum phase transition exists, but the high 
density phase actually shows no long range order in a strict 
sense. It exhibits a powerlike decay of the correlation func- 
tions both for the radial and Goldstone modes. We may 
still call this phase an "ordered phase" in a somewhat 
weaker notion: The rcnormalized order parameter pa{k) 
does not vanish, implying the distinction between Gold- 
stone and radial modes and several other features charac- 
teristic for an ordered phase. Also the order parameter po 
vanishes only asymptotically for fc ^ 0. For a system with 
a characteristic infrared cutoff kph ^ one can effectively 
observe order. A similar behavior has been found [l3| for 
classical phase transitions, e.g. the Kosterlitz-Thouless Q 
phase transition. 

For d = 1 one expects for the ordered phase a behavior 
similar to a Tomonaga-Luttinger liquid [l5| with dynamical 
exponent 2 = 1 and a correlation function 

(^*(r,r)^(0,0) -(vV2+r=^)-3. (83) 

The relativistic form of the propagator suggests that the 
"relativistic kinetic term" involving two 97--derivatives 
should not be neglected for low dimensions. We therefore 
will enlarge our truncation and include the coupling V in 
sect. I Vim This modifies the qualitative characteristics for 
the flow in the Goldstone regime for d = 1,2. For d = 2 
we will find that both po and po settle to constant values 
as A; ^ 0. The fixed point (C) disappears - it is an artefact 
of a too simple truncation. For d = 1 we find indeed a 
relativistic correlation function ([55)) with z = 1. The flow 
shows again a (shifted) fixed point (C), constant w and po 
and A ~ fc^ , po ~ k^- 

The qualitative new features induced by the coupling 
V limit the direct use of fixed point (C) in the simplest 
truncation (which neglects V). Nevertheless, the proper- 
ties of the flow equation (j72p remain interesting in several 
aspects. One concerns the "initial flow" before a substan- 
tial relativistic kinetic term ^ V has been generated. We 
discuss a few details of fixed point (C) for the system (17^ 
in appendix D. 

Let us finally briefly explore the behavior of eq. ([7^ for 
large w - details can be found in appendix F. For d > 3 
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one finds that the flow of po and A stops as w ^ k~'^ grows 
to large values for /c — > 0. Also the anomalous dimensions 
rj and rjs vanish. For d < 3 the flow of the combination 
Xw^^^ is attracted towards a partial fixed point. Again, the 
asymptotic behavior behavior for k —^ is characterized 
by constant po and A , rj 0. However, one now finds 
asymptotically vanishing X ^ S ^ k~^^ , rjs = 2{d — 
3). For d < 2 an initially very large value of w decreases, 
consistent with the attractor property of fixed point (C). 



VIII. CROSSOVER TO RELATIVISTIC MODELS 
FOR LOW DIMENSIONS 

For the Goldstone regime in d = 1 and ci = 2 an im- 
portant qualitative shortcoming of our simplest truncation 
becomes visible if we include the term with two time deriva- 
tives in an extended truncation 



Vv = -V 



b*di 



(84) 



A nonvanishing coupling V will always be generated by the 
flow of Ffe in the SSB regime, even if one starts with V = 
in the "classical action" at the microscopic scale A. This 
contrasts with the symmetric regime, relevant for the dis- 
ordered phase, where an initially vanishing V remains zero 
during the flow. For d — i the additional coupling V in- 
duces quantitative changes, but for small coupling the qual- 
itative changes in "overall thermodynamic quantities" , like 
density, pressure, order parameter and phase diagram, are 
moderate since the modifications of the infrared running 
only concern logarithms. Still, for more detailed features, 
like occupation numbers for small momenta, the coupling 
V is dominant. For the ordered phase in c? = 1,2, however, 
the relativistic dynamic term" (j84p will dominate over the 
term linear in dr and radically modify basic aspects of the 
macroscopic properties. In the Goldstone regime the cou- 
pling S vanishes for fc ^ such that the flow of the effective 
action is attracted to a (partial) fixed point with enhanced 
"relativistic" SO(d + 1) symmetry. This approximate rel- 
ativistic symmetry qualitatively changes the properties of 
fixed point (C). For d = 1 there will be a line of fixed 
points with different po ■ while the bare order parameter po 
vanishes ~ k"^ . For d = 2 the fixed point (C) disappears. 
The flow for /c — > will yield w — )■ oo and both po and po 
settle at constant values, with 77 = 0. 

We emphasize that the enhanced SO{d + 1) symmetry 
concerns only the leading dynamic and gradient terms for 
the Goldstone mode. It is not expected to become a sym- 
metry of the full effective action since the Lorentz sym- 
metry is not compatible with Galilei symmetry for T = 0. 
For example, an SO{d + 1) violating term with two time 
derivatives for the radial mode is possible, cf. app. G. 

For an initially vanishing or very small V a nonzero value 
is generated by the flow equation (A = XS^^k'^^^) 



dtV = 



av 



5-2 



(85) 



5vd 
d 



d + 2 



)a^(i + |) {l + wr^/^. 



(Details of the computation of the flow equation for V can 
be found in appendix C.) The relative importance of the 
kinetic terms linear or quadratic in dr can be measured by 
the ratio 

^ (86) 



kVV 

As long as s remains larger than one one may guess that the 
effects of S could remain important. Indeed, a naive scaling 
criterion for equal importance of the terms ~ S' or is 
given by Vqo w S with Sqo « k^ such that Vk"^ w S'^. We 
will argue, however, that for the Goldstone boson physics 
the relevant scale is \/2Apo rather than k. The effects of the 
couphng V therefore become dominant for V » S^/ (2Apo) 
or S* < y/w. 

For s ^ the effective action shows an enhanced 
SO{d + 1) symmetry, where t' — t/ \fv acts like an addi- 
tional space coordinate. From eq. ([55)) it is clear that the 
evolution of V essentially stops for fc ^ if S' decreases 
faster than k (and ay remains bounded). This will be the 
case for r\s < —1, but a weaker condition will be sufficient 
for an effective stop of the running of V. Indeed, from app. 
C we get the fiow equation for s 



dts = -(1 + Tjs)s + ^Av{s,w,X)s^ 



(87) 



where 



lim Ay — av 



hm A 



vs 



Xw 



One concludes that s is driven to zero if ?7s < —1- This 
presumably happens for d = I and d — 2. In this case 
the trajectories corresponding to an enhanced SO{d + 1)- 
symmetry are attractive - the long distance physics be- 
comes effectively relativistic. For 775 > —1 large values 
of s decrease and small values increase, suggesting a par- 
tial fixed point s^{X,w). If this occurs for large s we 
find ^ A^^/^i/;^/^. The relevant question for omitting 
the linear dynamic term ~ 5* in the Goldstone regime is 
s/y/w <C 1. This condition is reached for fc if V and 
S/ X go to constants, while S goes to zero. Constant values 
of V and S/X are suggested also on physical grounds since 
these quantities correspond to thermodynamic observables, 
see eq. (|G.27p in app. G. 

In this context we note that 5 = is always a (partial) 
fixed point, due to an enhanced discrete symmetry r — > — r 
(while keeping (j) fixed). (This additional discrete symme- 
try is preserved by our cutoff i?fe(16), even though this cut- 
off does not respect the SO{d + 1) symmetry - see app. C 
for a discussion on this issue.) For rjs < —1 the fixed point 
at s = is IR- attractive, while for 775 > —1 it becomes re- 
pulsive. For d = 3, where i]g > —1, the flow therefore ends 
for fc — !■ with nonzero s, corresponding to a violation of 
SO{d + 1) symmetry in the radial sector. For d = 3 one 
expects that V stops running for fc — > due to ui — s- 00. 
For large w one finds in eq. ([85]) ay = {X/2poy^^k^/S 
and dtV ^ S ^ (Infco/fc)^^. For A/5 const we therefore 
have logarithmic behavior 



A- S* 



1 

ln(fco/fc) 



(89) 
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and we note the difference as compared to tlie simplest 
truncation (|F.7p . where A decreases with the square of 
the inverse logarithm. This implies that s diverges ^ 
k~^\n{kQ/k) such that for d = 3 the large s regime ap- 
plies. 

Also for d — 1,2 the running of V stops, this time due 
to SO{d + 1) symmetry. The running of po and A within 
the Goldstone regime in the relativistic models has been 
intensively studied by non-perturbative flow equations [13] , 
[To| . For w ^ oo the running of po stops. On the other 
hand, the fluctuations of the Goldstone modes produce a 
fixed point for the dimensionless coupling Xk'^~^ for all d < 
3. One infers for the effective momentum dependence of the 
quartic coupling 



(90) 



Comparison with the simplest truncation (jF.14p shows that 
rjs has to be replaced by c? — 3 instead of 2(d ~ 3). This 
underlines again the crucial importance of the relativistic 
kinetic term for the long distance physics in all dimension 
d<3. 

The summary of the situation in the Goldstone regime 
is rather simple. For all d the asymptotic value for V 
reaches a constant as fc — > 0. For d = 3 also S becomes al- 
most constant (it vanishes only logarithmically), whereas 
for d = 1 and d — 2 the flow rapidly approaches an en- 
hanced SO{d -f l)-symmetry due to S vanishing with a 
power law S ~ k~'^'^ , S/k 0. The value of the renormal- 
ized order parameter po approaches a constant for d > 1. 
For d > 1 the anomalous dimension 77 vanishes for w —^ 00 
and also po become constant. The renormalized quartic 
coupling shows a scaling behavior according to its canoni- 
cal dimension in the relativistic model. 

For d = 1 and d = 2 the consequences of the "relativis- 
tic asymptotics" are immediate - the Goldstone regime is 
described by the classical 0(2)-model in 1 dimensions. 
With t' = T l\jv , gg = qoVV^ the correlation function for 
large distances in space and time (or small momenta go) 
obey (G = (Gil + G22)/2) 



G~(f2 + r'2)V 



G^(g^+g^^)-^+''/^ 

G^(f2+r'2)^. (91) 



(We recall that G is dominated by the Goldstone contribu- 
tion.) One may generalize the concept of dynamical crit- 
ical exponent also for situations without a finite correla- 
tion length. For d — 1,2 the effective dynamical critical 
exponent takes the "relativistic value" z = 1. For d = 2 
the Goldstone regime is described by the three-dimensional 
classical model. It is well known that po and po settle to 
constants, with 77 (fc — ^ 0) — 0. 

At this point we can already extend our discussion to an 
arbitrary number M of complex fields. The potential u{p), 
the gradient term and the relativistic dynamical term ~ V 
all obey an extended 0(2M)-symmetry. For our trunca- 
tion, the asymptotic behavior for the flow equations in the 
SSB regime is therefore well known for d — 1 and d — 2. 
Since S vanishes asymptotically, and S is the only term 



in our truncation that violates the 0(2M) symmetry, the 
asymptotic behavior of the flow is given by the classical 
0(2M)-models in d + 1 dimensions. (A more general dis- 
cussion of M-componcnt models will be given in sect. 1X1 ) 
In particular, for d — 2 one finds a simple description of 
order for arbitrary M in terms of the three-dimensional 
classical 0(2M) models. 

For d = 1 the two dimensional classical model applies. 
By virtue of the Mermin- Wagner theorem we know that 
no long range order exists with a spontaneously broken 
continuous symmetry. Since any po 7^ would lead to 
spontaneous breaking of the J7(l)-symmetry we can con- 
clude po{k — > 0) = 0. The way how the Mermin- Wagner 
theorem is realized depends on the number of components 
AI Q. For M > 1 both po{k) and po{k) reach zero at some 
positive value ksR. For kpt < ksR no order persists, while 
for kpfi > ksR the system behaves effectively as in the pres- 
ence of order. Typically, ordered domains exist with size 
Ld ^ ^sfl- Since the running of po is only logarithmic the 
scale ksR can be exponentially small. For an experimental 
probe with size L one has kph > so that for practical 
applications an "ordered phase" will persist. The typical 
size of ordered domains is then larger than the size of the 
system. (This issue has been discussed in detail for classical 
antiferromagnetism in two dimensions [2^.1 

For M — 1, in contrast, po reaches a constant value for 
fc ^ 0. Only the bare order parameter vanishes due to 
a nonvanishing anomalous dimension, po ^ fc**, such that 
order does not exist in a strict sense. In the correspond- 
ing classical model this situation describes the "low tem- 
perature phase" related to the Kosterlitz-Thouless phase 
transition. For practical purposes this phase behaves like 
an ordered phase, with powerlike decay of the correlation 
function G (PT|) due to the existence of a Goldstone boson. 
This is also the characteristic behavior of a Tomonaga- 
Luttinger liquid. It is well known from the classical 0(2) 
model in two dimensions that the low temperature phase 
is characterized by a line of fixed points which may be 
labelled by po = Po(k — > 0). The anomalous dimension 
depends on po [§, [lj| 



1 



V 



(92) 



It seems plausible that po depends on the effective chemical 
potential a such that we predict an anomalous dimension 
depending on cr. 

It is remarkable that the main qualitative features for 
d = 1 and fc — > 0, namely a nonzero po, vanishing po, and 
a positive anomalous dimension rj > 0, are already visi- 
ble from fixed point (G) in the simple truncation of sect. 

Not surprisingly, however, the quantitative accuracy 
for the anomalous dimensions is very poor if the coupling 
V is omitted. We may indeed address the properties of the 
Goldstone regime in the perspective of the properties of 
fixed points in presence of the coupling V. For d = 2 one 
has the well known Wilson-Fisher fixed point of the three 
dimensional classical model. It corresponds to 5* = 0. The 
question of how close trajectories approach the Wilson- 
Fisher fixed point depends on the microscopic parameters a 
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and Aa as well as on a possible microscopic coupling V{A). 
Quantum phase transitions with critical behavior different 
from eq. ([52|) can be associated with the Wilson-Fisher 
fixed point. In this case z = 1 and the critical exponents 
1/ and rj of the three-dimensional 0{2M) model apply. For 
d = 2 this type of phase transition presumably becomes 
relevant for large enough microscopic couplings V{A). For 
V{A) = 0, as considered in this paper, the quantum critical 
fixed point discussed in sect. V is relevant. For this quan- 
tum critical fixed point a vanishing relativistic coupling 
V — is stable with respect to the flow. In our truncation 
we infer from eq. (|85p that for w — one has ay = and 
therefore dtV = while there is anyhow no contribution to 
dtV in the disordered phase. At the quantum critical point 
the dimensionless combination Vk^ therefore corresponds 
to an irrelevant coupling. 

In order to judge the relevative importance of the 
Wilson-Fisher (WF) and the quantum critical (QC) fixed 
points for arbitrary microscopic couplings V^(A) one should 
consider the critical hypersurface on which both fixed 
points lie. (Note that po{A) varies on this hypersurface, 
with po{A) = for QC and po(A) > for WF. We use a 
common name (QC) for fixed point (A) {d > 2) or fixed 
point {B) {d < 2)). The first question concerns the sta- 
bility of WF with respect to the coupling S. Taking into 
account the scaling dimensions at WF one finds that WF is 
stable for 775 < —1 and unstable for 77s > —1. Here r]s has 
to be evaluated for WF, which we have not done so far. For 
ijs > —1 one would observe a crossover from WF to QC 
on the critical hyperface. In contrast, for rjs < —1 both 
WF and QC are stable on the critical hypersurface. The 
topology of the flow would then imply the existence of a 
new fixed point with finite nonzero value of S{A)/ ^yv{A). 

For d = 1 (and M = 1) the role of the Wilson-Fisher 
fixed point is replaced by the Kosterlitz-Thouless fixed 
point for the two dimensional classical 0(2) model. A 
key new ingredient is the existence of a whole line of fixed 
points for S = 0. They can be parameterized by the renor- 
malized order parameter po (corresponding to k in Ref. [14] 
at fc = 0. These fixed points govern the Goldstone regime 
of our model with V{A) = 0. Thus the IR attractive fixed 
point (C) in the truncation with y = transforms into one 
of the fixed points on the critical line. Now w is no longer 
an irrelevant coupling - it can be used to parameterize the 
line of fixed points instead of po. (Indeed, w = 2(A/fc^)/9o 
and (A/fc^) approaches a fixed point value depending on pQ 
0.) It seems natural that po depends on a. On the other 
hand, poi^k = 0) cannot take arbitrary small values, corre- 
sponding to the jump in the renormalized superfluid den- 
sity of the Kosterlitz-Thouless transition [23] . This raises 
interesting questions of how the chemical potential a is 
mapped into an allowed range of po or w. It is likely that 
the answer is linked to the "initial flow" for small V with 
a possible influence of an approximate fixed point of type 
(C) for which is a small perturbation. 

We observe that on the line of fixed points the anomalous 
dimension 77 depends on pQ. The maximal value rj — 1/4 is 
reached at the Kosterlitz-Thouless transition. Our trunca- 
tion therefore predicts ?7 < 1/4 for eq. (pij) . and we infer 



for d = 1 , M = 1 

G~ (r2 + r'2)-^. (93) 

However, it is known for the Tomonaga-Luttinger liquid 
that the maximal value for 77 is 1/2. This brings us to the 
question of reliability of a truncation that only includes a 
potential, second order gradient terms and first and second 
order time derivatives with SO{d-\- 1) symmetry for S = 0. 
For V > such a truncation is not compatible with Galilei 
invariance for T = and p ^ 0. We conclude that the 
description in terms of a d-l- 1-dimensional 0{2) model be- 
comes problematic for w < 1. On the other hand we have 
presented strong arguments for its validity for w ^ 1. We 
conclude that the identification of the Tomonaga-Luttinger 
liquid with the Kosterlitz-Thouless phase of the two di- 
mensional model may only be valid for large w. For w ~ 1 
certain terms neglected in our tuncation and breaking of 
the SO{d -\- 1) symmetry can play a role. The range of 77 
near 1/4 in the two-dimensional 0{2) model corresponds to 
w 1 - this is the region of the phase transition. We con- 
clude that our approximation is expected to hold as long 
as rj is small enough (large w), but may break down for 77 
near 1/4, such that the upper bound rj < 1/4 is not robust. 

We finally comment on the case d = 3 where the lin- 
ear T-derivative remains important, as far as naive scaling 
is concerned. However, in the Goldstone regime the lin- 
ear T derivative involves the radial mode ipi, J (p*dr'p ~ 
J (pidr^2- If we study the long range physics the radial 
mode effectively decouples for distances larger than the cor- 
relation length and we should question the criterion s < 1 
that we have used for an estimate of the importance of the 
term ~ V. This criterion is valid for physics involving the 
radial mode, i.e. kph We will discuss in the next 

section that for kph ^ Cij"'^ the relevant dimensionless ratio 
is not given by s but rather by S/y/2XpoV. This can be 
seen by "integrating out" the radial mode and investigat- 
ing the effect of 5 in a nonlinear model for the Goldstone 
boson. 



IX. GOLDSTONE REGIME AND NON-LINEAR 
cr-MODELS 

In the Goldstone regime the infiuence of the radial mode 
is subdominant. We may therefore aim for a description 
only in terms of Goldstone bosons. This leads to the non- 
linear cr-niodels. The effective action in the nonlinear for- 
mulation can directly be computed from our results. One 
solves the field equation for the radial fluctuations as a 
functional of the Goldstone fluctuations and reinserts this 
solution into the effective action. For this purpose we pa- 
rameterize 

(t,{x)^{cf>o + R{x))e''(-^ (94) 

and eliminate the radial mode R{x) as a functional of the 
periodic phase variable 9{x). 
Within our truncation 

I C, C^u{(j)*(t)) + S(j)*dr(t>-(j)*{A + Vd^)(t> (95) 
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we find (for constant S, V and removing total derivatives 
and constants) 

C = Co + Cr, 
Co = 4>1{V0V0 + Vdr0dre}, 
Cr = u{<l>l + 2(j,oR + R^)-u{<l>l) 

+ (2(/)oi? + R'^){V0V9 + Vdr0dr0 + iSdr0} 
+VRVR + VdrRdrR. (96) 

The field equation SCr/SR = can be solved iteratively 
by expanding in {V0V0 + Vdr0dr0 + iSdr0}. In lowest 
order one obtains, with A = ^"(^o), 

R = -M'^>'<l}l-^-V^d^y^{W0W0 + Vdr0dr0 + iSdr0}, 

(97) 

and reinserting this solution yields (po = 0o) 

Cr = -po{iSdr0 - Vdr0dr0 - V 0V 0} (98) 
(2A/90 - A - Vd^)-^{iSdr0 - Vdr0dr0 - ^0W0}. 
Expanding in powers of 



Cr, — Cr.2 + Cint, 

Cr,i - poS''dr0{2\po 



Vdl)-^dr0 (99) 



and combining with Co we obtain the effective inverse prop- 
agator for the Goldstonc boson (^ 4>o0) 



9 



(100) 



On the other hand, Cint contains interaction terms. These 
are purely derivative interactions, as appropriate for Gold- 
stone bosons. Higher orders in the iterative expansion will 
only yield further interaction terms. 

The inverse Goldstone propagator shows no term linear 
in qo , even for V = , S = \. This is a simple consequence 
of the observation that a linear r-derivative for a single 
real field is always a total derivative and can therefore be 
eliminated from the field equations. In lowest order in a 
derivative expansion (for small momenta) one has 



Gg — 



V 



2Xpo 



(101) 



such that both V and S contribute to an effectively rel- 
ativistic kinetic term. (In the limit of large go i ^ 
2Xpo + q the correction from Cr results in an effective 
masslike term S'^/V, i.e. Gg = [q"^ + V(g + S"^ IV\-^ :) As 
mentioned at the end of the preceeding section the relevant 
ratio for the importance of the linear dynamic term ~ 5 is 
given for the Goldstonc regime by Sj \J2\poV ■ 

We conclude that in the Goldstone regime the renormal- 
ized propagator always takes a relativistic form 
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The quantity 



v'=\V 



2Apo 



(102) 



(103) 



can be associated with the microscopic sound velocity 
of the Goldstone mode. Comparing with eq. (|G.23p 
in appendix G we see that it equals the macroscopic 
sound velocity, Cg in accordance with Ref. [3]. For 
po const, S/X const, S ^ 0, the second term 
in eq. p03p can be neglected. This yields a simple relation 
between the sound velocity and the coupling V = V{k 0) 



The bare propagator reads correspondingly 



Ga 



[q' + iqoMY^'-^^ 



(104) 



(105) 



where rj depends on momenta according to the replacement 
q^ +q^/v'^ . For d> 1 the anomalous dimension van- 
ishes asymptotically, ri{k ^ 0) = 0. The generalization to 
M component fields is straightforward and will briefly be 
addressed in the next section. We emphasize that for the 
nonlinear models the propagator shows no qualitative dif- 
ference between a linear dynamic term ~ 5' or a relativistic 
dynamic term V . Nevertheless, the form of the dynamic 
term becomes important for the linear model, in particular 
for the flow of the order parameter po{k). 



X. M-COMPONENT MODELS 

For models with a relativistic kinetic term the critical ex- 
ponents depend sensitively on the number of components 
of the field. For N real components and scalar models 
with 0{N) symmetry the universality classes can be fully 
characterized by N and d. We may investigate this ques- 
tion also for the nonrelativistic kinetic term and investigate 
models with M complex scalar fields , m — \ . . . M. 
In terms of renormalized fields we consider models with 
[/(M)-symmetry and neglect first the coupling V 

Tfc = / { Y^^S^rndrCl^"^ " CA<^™) + u{p)} . (106) 

Due to the U{M) symmetry the potential can only depend 
on the invariant p = '^4'm4'm- Since p is invariant with 

m 

respect to the symmetry 0(2M) acting on the 2M real 
components of the scalar field, the potential exhibits an 
enhanced 0{2M) symmetry. This enhanced symmetry is 
shared by the term involving spatial gradients but violated 
by the linear r-derivative. 

Within our truncation of pointlike interactions and with 
a linear r-derivative, we find a very simple result: the fiow 
equations do not depend on AI, except for the flow of a field 
independent constant in u. Neglecting the momentum de- 
pendence of interactions all our discussion therefore holds 
without modifications for the M-component case as well. 
As long as the propagator is specified by eq. (|106p a de- 
pendence of the critical exponents on M can therefore only 
arise from the momentum dependence of the interactions. 
This situation is quite different in presence of a relativistic 
kinetic term involving two r-derivatives. 
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In order to show this result we first note that for a con- 
stant background field the inverse propagator is block di- 
agonal, G~^= diag{Gi^ , Gj"^ . . . ) where are 2x2 
matrices corresponding to a decomposition of the complex 
field (jjm into real and imaginary parts. Without loss of 
generality we may take a real constant background field 
in the m — 1-direction. Then G^^ is given by eq. ((20|) 
whereas for m > 1 one has 

G-u - i V ' , ] • (107) 

\ Sqo , q +u' J 

The resulting correction to dtu from the (M — 1) additional 
components does not depend on p 

A field independent additive constant does not influence 
our discussion and we recover the same flow equations for 
w and A as for M — 1. 

We next turn to the computation of ij and 775. We define 
these quantities by the flow of the q ^ and go dependence 
of the inverse propagator for the m = 1 component, ac- 
cording to eqs. ([551) . We can use eqs. ([HO)) . ([^ with 
indices a,b . . . , running now from 1 to 2Af , e.g. the m = 2 
component corresponds to a = 3,4 etc.. In presence of a 
background field in the a — 1 direction the interaction term 
in the effective action for the fields (j)b,b 1 has a discrete 
symmetry — s- —0b. This implies that for nonvanishing 
cubic couplings "facd the values of all indices must appear 
in pairs, except for the value one. For the anomalous di- 
mension r] the indices a and b in eq. pip take the value 
two. In consequence, the index pair (d,e) can only be (1,2) 
or (2,1), and similar for the index pair (/, c). Therefore no 
indices c,d,e,f with values three or higher appear and the 
additional {M — 1) complex fields do not influence rj. For 
775 the situation is similar if one also uses the fact that G 
is block diagonal, i.e. Gi^b>2 = G2,b>2 = 0. (The situa- 
tion can be easily depicted in terms of one-loop Feynman 
graphs with two external legs.) This closes our argument 
that the flow equations are independent of the number of 
components M in the approximation of pointlike vertices 
and for V = 0. 

We conclude that the symmetric phase and the quan- 
tum phase transition are not influenced by the presence 
of additional components of the field. At the quantum 
critical point one finds for all M the critical exponents 
1/ = 1/2 , 77 = , z = 2. In the ordered phase the number 
of components will influence the flow as soon as a sizeable 
strength V for the second r-derivative is generated. In par- 
ticular, the flow for d = 1 and d = 2 will depend on M. 
The flow for fc — > is expected to reproduce the well known 
behavior of the Goldstone boson physics in classical statis- 
tical systems H, [l^] ■ Typically, the asymptotic behavior 
can be described by nonlinear a-models on a manifold given 
by the coset space C/(Af)/[/(M - 1). For d-h 1 = 2 the be- 
havior of the abelian model for M = 1 differs substantially 
from the non- abelian cr-models for M > 1. 



In the microscopic action only the coupling S violates 
the 0(2M) symmetry. In analogy, we expect that for 
small 0(2M) violating couplings, i.e. small S we can 
choose a description in terms of 0(2M) /0(2Af — l)-models. 
This will contain symmetry breaking effects since the 
0{2M) symmetry is exact only for S* = 0. Such ef- 
fects are proportional to the dimcnsionless combination 
S{k)/{y^2X{k)po{k)V{k)). Since for d = 1,2 the flow is 
attracted towards S = (if t]^ < 0), one may suspect that 
the enhanced symmetry for 5* = 0, i.e. the space-time ro- 
tations 50(^-1- 1) and the internal rotations S'0(2M), are 
approximately realized in the Goldstone regime. The de- 
gree of violation of these symmetries depends on the char- 
acteristic momentum scale of the Green's functions and on 
w~^. As before, a nonzero "external momentum" q may 
be associated with fc^^ = q"^ > 0. Since 5(fc) vanishes 
only asymptotically for fc — ^ the symmetry breaking of 
50(^+1) and 50(2M) due to S{k^h) > will always be 
present. Furthermore, the vanishing of S occurs for large 
w, while S remains important for w < 1. 

The order parameter po oi' Po is a quantity that in- 
volves the limit of zero momentum (or kph « L^^ with 
L the macroscopic size of the probe). We may there- 
fore take the limit fc ^ , g ^ — > for the issue of 
spontaneous symmetry breaking. For small momenta we 
may consider a derivative expansion of the effective ac- 
tion (for fc — > 0). Let us consider up to two derivatives. 
For 1 < d < 2 both w^^ and the ratio S{k)/k vanish for 
fc ^ and s can be neglected. The potential and the rela- 
tivistic two-derivative term —ip*{VdT'^ + A)(p respect both 
SO{d + 1) and 0(2A/). A term Yp{VdT^ + A)p vio- 
lates SO{d + 1) only iiV j^V, while it preserves 0(2M). 
On the level of two derivatives an 0{2M) violation could 
only arise from a term ^^{if\ipidr'P2dT^2 + V2V2dT^idT^i 
~(p2^idT^idT^2 — Vi'P2dTiP2(^TVi) for M = 2, or suitable 
generalizations for M > 2. (On the level of two derivatives 
the invariants preserving U{M) while violating 0{2M) are 
of the form | £„!... mjv'^mi • ■ i9(/),„„ p.) 

For d = 2 we may speculate that the SO{d + 1) and 
0(2M) violating terms involving two derivatives vanish for 
fc ^ 0. The behavior of the order parameter and the corre- 
lation function at small momentum are then well approxi- 
mated by the classical 0{2M) model in three dimensions. 
This issue depends on whether the Lorentz symmetry vi- 
olating operators are irrelevant for the trajectories with 
w 00, or not. Even in presence of 0(2M) violating 
derivative terms one would still find 2M — 1 gapless ex- 
citations, due to the accidental 0(2A'/) symmetry of the 
potential. 

For d — 1 the issue may be more complicated. Within 
the relativistic model w is either constant (for M = 1) 
or runs logarithmically towards zero (M > 1). We may 
therefore encounter a region where w < 1 such that the 
SO{d + 1) violating terms can no longer be neglected. For 
M > 1 this region becomes always relevant for fc — > 0. 
For M = 1 we expect that SO{d + l)-symmetry becomes 
a good approximation for the fixed points with small 77 
(large w), whereas it is questionable for the larger values 
of 77 corresponding to w ^ I. 



18 



XI. CONCLUSIONS AND OUTLOOK 

The functional renormalization group yields a unified 
picture for quantum phase transitions of bosons. It covers 
both the ordered and the disordered phase, including the 
rather complicated long range dynamics of the second. We 
describe within the same simple truncation an arbitrary 
number of space-dimensions d and an arbitrary number 
of components M (for M complex fields with symmetry 
U{M)). In a general setting, the properties of the phase 
transition depend on the parameter sa which characterizes 
the relative strength of the kinetic terms involving one or 
two time derivatives. More precisely, for a general micro- 
scopic (or classical) action ([T]) the dimensionless ratio sa is 
given by 



SA 



2Mb 
V 



SiA) 



(109) 



with A the ultraviolet cutoff. 

The characteristic features are described by two limits. 
The point sa = describes models with a "relativistic 
kinetic" term involving two time derivatives and an en- 
hanced space-time symmetry 50(^-1- 1). For those models 
the quantum phase transition in d dimensions is strictly 
equivalent to the classical phase transition in c? -t- 1 dimen- 
sions. The universal critical properties correspond to the 
well studied 0(2Af )-models in 1 dimensions. 

For the other limit sa ^ oo (i.e. = 0) Galilei symme- 
try characterizes the zero temperature physics. Now the 
phase transition is influenced by the non-relativistic quan- 
tum critical point. This is the case we have mainly studied 
in this paper. For intermediate values of sa one expects 
near the phase transition a crossover from the "classical" 
or "relativistic" critical point to the non-relativistic quan- 
tum critical point. No phase transition is expected as sa is 
varied. However, the critical exponents and amplitudes will 
depend on sa. They are given by the universality class of 
the d+1 dimensional relativistic 0(2M) models in the limit 
of small Sa - for example by the well known Wilson-Fisher 
critical point for d = 2. As sa increases towards infinity 
the critical exponents smoothly change towards the mean 
field values for the non-relativistic quantum critical point. 
We note that for sa ^ the critical hypersurface for the 
phase transition typically occurs for nonzero density. 

The critical behavior for the nonrelativistic "quantum 
critical models" (i.e. sa — oo) can be characterized by two 
regimes, depending on l/D. Here I is a characteristic length 
scale of the system or experiment, and D is proportional to 
the average distance between (quasi-) particles, D ^ n~^^'^. 
The "dilute regime" applies to the disordered phase (where 
D oo) and to the "linear regime" in the ordered phase, 
as shown in fig. [1] In the dilute regime simple scaling laws 
apply and the critical behavior is described by mean field 
critical exponents for all d and M. The simplicity of its 
properties finds a simple explanation, since the disordered 
phase describes the vacuum with zero particles, and for the 
linear regime in the ordered phase the small density gives 
only subleading corrections. In more technical terms, this 
is the regime where fc^^ = l^^ ^ 2Xpo or w <C 1. 



For I ^ D the particle density matters. For this "dense 
regime" the particle density n sets a new relevant mo- 
mentum scale ~ or energy scale {2MbD^)~^ . The 
long distance physics is described by the interacting Gold- 
stone bosons which arise from the spontaneous breaking of 
[/(M)-symmetry. We therefore often call the dense regime 
the "Goldstone regime". Since Goldstone bosons must be 
massless (or "gapless" ) the correlation functions decay with 
an inverse power of distance in space or time. 

The physics of interacting Goldstone bosons crucially de- 
pends on the number of space dimensions d. For d = 3 we 
find that the "Goldstone dynamics" is influenced both by 
short and long wave length fluctuations. The running of 
dimensionless couplings is logarithmic. For small or mod- 
erate interaction strength Aa mean field theory remains a 
good guide. Quantum corrections induce quantitative cor- 
rections but do not change the qualitative behavior, except 
for the extreme infrared. We do not address in this paper 
the possibility that large Aa may lead to new phenomena, 
as the destruction of the condensate by a too large repul- 
sion between the bosons. 

On the other hand, for d = 2 and d — 1 the Goldstone dy- 
namics is infrared dominated, leading to qualitatively new 
features induced by quantum fluctuations. One expects 
strong deviations from mean field theory. We find that the 
renormalization flow describes a crossover to models with 
a relativistic kinetic term. Even for models with sa — oo 
(i.e. = 0) the value of s{k) decreases fast for momentum 
scales k A. The relativistic term ^ V will be induced by 
the fluctuations and dominates for k 0. The Goldstone 
regime for d = 1, 2 is therefore characterized by s = 0, both 
for Sa — *■ oo and for sa = 0. In other words, the flow is 
attracted fast towards the (partial) fixed point s = 0. (The 
flow in the Goldstone regime differs from the flow on the 
critical hypersurface where s presumably diverges.) 

This implies for d = 1 , 2 a close link between the Gold- 
stone regime of the quantum model in d dimensions and 
the corresponding Goldstone regime of the classical model 
in d + 1 dimensions. The classical Goldstone regime is 
well studied. For d — 2 one expects a description of the 
Goldstone regime by the classical three-dimensional U{M)- 
model. In this case the realization of order is straightfor- 
ward, with both po and po approaching constants for A: — > 
and rj tending to zero. The correlation function for large 
separations in space or time shows a simple powerlike be- 
havior 
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Here we have restored dimensions, with sound velocity Cg = 
v/{2Mb). 

For d = 1 the "ordered phase " actually shows no long 
range order in a strict sense. The influence of the Goldstone 
bosons is so strong that they destroy any long range order. 
This agrees with the Mermin- Wagner theorem for the clas- 
sical two-dimensional models with continuous global sym- 
metry, as applicable here. The way how this theorem is 
realized depends crucially on M . For M = 1 the non-linear 
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cr-model is abelian. In the ordered phase the relevant ex- 
citations include vortices. Indeed, the Kosterlitz-Thouless 
phase transition finds within the functional renormaliza- 
tion group a simple description in terms of a linear 0(2) 
model in two dimensions ,14]. As a result, the renormalized 
order parameter po reaches a constant for fc — + 0, while the 
bare order parameter po vanishes ~ due to a nonvanish- 
ing anomalous dimension. The correlation function decays 



(Ill) 



These findings carry over to the non-relativistic model for 
d = I, the Tomonaga-Luttinger liquid. In our approach the 
key features of the Tomonaga-Luttinger liquid result essen- 
tially from the nonzero value of the renormalized order pa- 
rameter Po- This model simply describes the Goldstone bo- 
son associated to the "effective U{1) symmetry breaking". 
The anomalous dimension depends on the density of quasi- 
particles, n — po- For large k — \fV n — n/v — n/ {2MbCs) 
one finds from eq. ([^ 



MbCs 
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(112) 



After restoring dimensions the sound velocity Cs is related 
to the couphng V by eq. (|G.27p 

(113) 
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For M > 1, as for example for the ferromagnetic to para- 
magnetic transition with M = 3, the nonlinear a-models 
are asymptotically free and induce a nonperturbative scale 
ksR- Within a description by linear U{M) models one finds 
that po{k) and po{k) both vanish for k = kg^, while be- 
ing finite for k > ksR- For k < ksR the flow follows the 
symmetric regime with a minimum of the potential at the 
origin. The situation with nonzero po{k) describes local 
order, with a maximal size of domains up to kg^. The 
order in larger domains is destroyed by the Goldstone bo- 
son fluctuations. Since ksR can be exponentially small for 
large enough density (the running of poik) is only logarith- 
mic) , there is always a critical density ric beyond which the 
size of ordered domains fc^^ grows beyond the macroscopic 
size of the experimental probe L. Thus for n > Uc one ob- 
serves effectively spontaneous symmetry breaking, despite 
the Mermin- Wagner theorem. 

The dynamical behavior of the quantum critical models 
can be characterized by a dynamical critical exponent z. In 
general terms, it describes how a characteristic time scale 
f scales with a characteristic length Z, i.e. f ~ P. We 
have concentrated on the non-relativistic quantum critical 
models (sa — > oo). We find for all d and M that the 
value of z depends on the density of (quasi)-particles. For 
the dilute regime one has the mean field scaling z = 2, 
while the dense regime shows the relativistic value z = 1. 
Inbetween, there is an effective crossover. As far as time 
scales are concerned the Goldstone regime (dense regime) 



applies for f > D/cs, with D the interparticle distance 
~ ^-i/rf and Cg the sound velocity. 

This paper has demonstrated that the flow of a small 
set couplings, namely pQ (or m?), X, S and V yields all 
qualitative features for the quantum phase transition for 
arbitrary d and M. The quantitative precision can be im- 
proved by extending the truncation. Straightforward steps 
include the incorporation of a nonvanishing V in the flow 
equations for the effective potential u and for rjs, as well 
as a differentiation between the effective renormalization 
constants for the radial and Goldstone modes via the in- 
clusion of a term Ydpdp. Together with 7 = "'•'^^(po) we 
expect that the set of couphngs (po. A, 7, S, V, Y) will yield 
already a very satisfactory quantitative accuracy. Exten- 
sions to include fermionic degrees of freedom are possible 
[25j . Furthermore, the effects of nonvanishing temperature 
can easily be incorporated in our framework - one sim- 
ply has to replace the go-integration by a Matsubara sum. 
Close to the critical temperature of the phase transition 
and away from the quantum critical point at T = the 
flow will experience an effective dimensional reduction to 
the classical c?-dimensional U{M) models. High accuracy 
for the functional flow equations has already been reached 
for the classical 0{N) models in arbitrary dimension. In- 
cluding the temperature effects the flow equations should 
provide a rather complete picture for the thermodynamics 
of bosonic quantum gases. 

In this paper we have concentrated on the structural 
aspects. A numerical solution of the flow equations for k 
will directly yield important thermodynamic quantities 
in the zero temperature limit. The density is given by 
n = Po and the condensate fraction by fie = A~^. The 
sound velocity obeys Cs = 1/{2MbW) for d = 1,2, with 
logarithmic corrections for d = 3 due to S ^ 1/ ln(fco/fc), 
cf. app. G, 
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The change of the condensate fraction with the density 
follows (cf. app. G) 
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with A the renormalized quartic coupling. 

It may be possible to measure all these quantities by 
investigating the Bose-Einstein condensate for ultracold 
bosonic atoms. With suitable traps one may prepare essen- 
tially homogeneous systems for d = 1, 2 or 3. For a quanti- 
tative computation one will further need to determine Aa 
in dependence on external parameters, as a homogeneous 
magnetic field. For this purpose one relates it to measur- 
able properties in the vacuum, like the scattering length for 
d = 3. The necessary computation for n = corresponds 
precisely to the phase transition discussed in sect. IVl and 
can be performed in the comparatively simple disordered 
phase. The prospects for experimental tests of some of the 
features described in this paper look promising. 



20 



Note added: 

After the first version of the paper interesting extensions 
and numerical results have been obtained in Ref. 12611. 
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The second derivative reads 



APPENDIX A: Flow equation for the effective 
potential 



We can write the flow at fixed renormalized field (j) 
dtu — 7]pu' (A.l) 

+ / {[k''l{e-ci')\e{e-<i')trG] 

k 



where v^^ = 4:n,V2 ^ — 87r,t;3 ^ — Stt^ and q — {q^Y^^- 
The integration over qo can be carried out easily 



2. 27r 



dqo 



q +u' + pu" 



(27r) ^2^2 ^ (~2 + y,)(^2 + + 2pu") 

I q^ +u' + pu" 



25 ^q^ + u'^q^ +u' + 2pu" 



(A.2) 



In the integrand we can use cf' = fc^ such that the q inte- 
gration is trivial and yields eq. ((22| . 



Differentiation of eq. ([22]) with respect to p yields 



dtu' = ri{u' + pu' 
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dS 
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■q 



d + 2 
I 



VWTu' ^fc2 +u' + 2pu" 
j pu"^ 3pu"^ + 2p2u"u(3) 



fc2. 



fc2 +u' + 2pu" 



(A.3) 



dtu" — r](2u" - 








d + 2 
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I 



VWTv' sjk'^ + m' + 2(m" 

\ m"2 + 2pv!'u'^^^ 
I k'^ +u' 

3m"2 _^ 10p^"u(3) + 2p2(u(3)2 + u"u(^^) 
k^ +u' + 2pu" 



(A.4) 
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/'2 (3u" + 2pm(3))2 



(fc2+M')2 (fc2 + m' + 2pu")2 

while the third derivative becomes already quite lengthy 



with 
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_2va ( _ d+2 
dS \ d + 2 J 
1 1 
VWTu' ^A:2 +u' + 2pu" 



R (A.5) 



fc2 + u' k + u' + 2pu" 



[k'^+u'f (fc2 +U')(fc2 +u' + 2pw") 

(fc2 + u' + 2pu"Y 



{k'^+u'f (fc2+u')2(fc2+M' + 2pM") 

23 

(fc2 +w' + 2pit")2 

Z4 



and 



(fc2 + u' + 2pu")3 
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+2p2 (3^(3)^(4) 5))^ 



Z4 



= (37/" + 2p«(3))(9^"2^^^^,,^(3) 
+6p2^(3)2^g^2^»^(4))^ 

= ^^pv!'\ Z2 = \pu!'^{ivl' + 2pu^''\ 
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pM"2(3u" + 2pu(3))2 ^ 



pu"(3w" + 2pu(3))3^ 



(A.7) 
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APPENDIX B: Computation of anomalous 
dimensions 



For the calculation of rj and dt In S we infer from eq. ([?T|) 



dtP22{q) = Piuyi' / dtRkiq') 

Jq' 

{iG')lliq')G22{q' + q) 

+(G2)22(g')Gii(g' + <z) 

+ {G^)l2{q')G,2{q' + q) 

+ {G^)2i{q')G2iiq' + q) + (q 



(B.l) 



-q)}, 



and 



(G2)ii(g')Gii((z' + g) 



+ {G^)22{q')G22{q' + q) 



+ (3 + ^) iG')uiq')G2M + q) 



— ] {G')2l{q')Gi2{q' + q) 



+iq^-q)}, 

dtPl2{q) - p{u")^A' f dtRkW) 



(B.2) 



{ (3 + ^^) {G')ll{q')[Gi2{q' + q) + G2i{q'-q)] 

+ {G^)22{q') [G2i{q' +q)+ G^q' - q)] 
+ {&)i2{q)G22{q +q) + (0^)21 (9)^22 (<z' - q) 

+ (3 + ^) [{G')2liq')GMq' + q) 

+ {&)l2{q')Gn{q' -q)]}, (B.3) 

5^^21(9) = p{u"fA' f dtRkiq') 



{ (3 + (G^)ii(g') [G2i(g' + g) + Gi2(g' - g)] 

+(G2)22(<z')[Gi2(g' + g) + G2i(g' - g)] 

+(G2)2i(g')G22(g' + g) + (G')i2(g)G22(g' - g) 
2pu(3) 



[(G2)i2(g')Gii(g' + g) 
+ (G^)2i(g')Gii(g'-g)]}. (B.4) 

The propagator matrix reads explicitely 



Sqo 



G = det-iA-i ( 9'+"' ' I , (B.5) 

I -5go , q^+u' + 2pu" ^ ' 



with 

G^ = det-^A"2 
/ {e+u'f-S^ql 



(B.6) 



25go(g2+w' + pw") 
25go(g2+M' + pu") , {if +u' + 2pu"f - S^ql) 



and 



det = + u')if +u' + 2pu") + S'^ql 



(B.7) 



For the computation of 9tP22 we employ the fact that only 
q'^<k^ contributes in the integrands (|B.4p and replace 
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{q'±qf 



z±. 



(B.8) 



Including terms up to second order in q we can expand 
in z± 

9*^22 (g,0) = AVo / / 0{k^ - q' ^)dt[A{e - q' 2)] 
2z+(fc2 + u' + pu") + z\ 
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dett 



with 



4z2(fc2+u'+pu")2 

det?. 



(fc2 + lt')2(fc2 + „' + 2pv!' + Z+) 

+ (fc2 + 7/ + 2pv!'f{k^ + li' + z+) 
+25'g^,2(fc2 + u' + pv!' - + (z+ ^ 



(B.9) 



(B.IO) 



and 



detfc = (fc2 + M')(fc2 + + 2pw") + S'2g^2_ (g-^^^ 



Here we use the fact that z± is effectively linear in q for 
(f 2 ^ since only momenta q' ^ « fc^ contribute for terms 
involving powers of 2+ or z- . Without loss of generality we 
choose (f = (g, 0, . . . ), g > and decompose q'^— q'^ + xt 
(for d> \). The combination of the ^-functions from dtRk 
and from z+ restricts the integration range for q ' to 



Q <xt <k 
max 



(B.12) 



{- - Xt-, \/k'^ - Xt - g} < g'l < \Jk'^ - 



Xt. 



In consequence, the term q'"^ — k^ in z+ is effectively of the 
order q. It will be convenient to split the q '-integration 
into two ranges 



/ : < < fc^ - 



\/k'^ - Xt — g < g'l < \/k^~~xt, 
-a^ <xt < fc^ 



k 



- \Jk'^ - Xt < gi < Vfc2 - 



(B.13) 
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For d = 1 there is no a;(-integration and the q[ integration 
covers the range k — q < q[ < k. Restricting the Xt and 
q'j-integrations to this range we write 

„ a - _ „ p{u"fA 



[(2 - 77)fc2 + r/gf + ijxt] det^ + ^2^+) 

+ ((Z^-g)| (B.14) 



with 



and 



Ai = detfc, 

A2 - -{k' + u' + pu") 



(B.15) 



z+=qf +a;t-fc2 + 2ggi+q2. (B.16) 

For d — 1 the Xf-integration and the factor 2vd-i are ab- 
sent. Using a = \/k'^ — xt and 



dq[ [(2 - 77)fc2 + ijxt + rjq[^] {AiZ+ + A^z^) 



= 2fc2g2yp~:^^^ (B.17) 
we obtain for d = 1 at the minimum p — poi ^^"(/Oo) = 



d 



2X^PoAk^ 



dt^P22\,^,^- ~" Zr I dq'odet-\ (B.18) 



For d > 1 we stih need to perform the Xj-integration and 
to include the region II. We employ 



K 

J dxtF[xt) = \q^F{e) (B.19) 

and observe that the integration region II does not con- 
tribute in order q^ . This yields for d > 1 



dt-Q^P22\q=o 



dq^ detfc 



(B.20) 



with 



We collect the identities (n, > 1 

1 



dqo detfc 



(n- 1)! 



(l-3-5...2n-3)^B"^, (B.23) 

00 00 

Jdq'oqo'det-^-+'^ = l- J dgWetr , 



with 



B = {k-" + u'){k-' + u' + 2pu") (B.24) 



such that 



(B.25) 



In terms of w and A we obtain our final result eq. 

For the computation of dtS we expand in linear oder in 

dtP2iiq^0,qo)^2p{u"fSqo 



dik'-q'')dt[Aik'~q")] 



9 ' 
det: 



2J2 i ^ , P'^ 



-2S'cl^' 1 



(3) 



} 



(B.26) 



We evaluate eq. (jB.26[) at the minimum v! = 0, tt" = 
The q '-integration is trivial and the q^- 
integration follows from (jB.23p . One finds 

[8 - 4w - 3^2 + (8 + w)po7/A] . (B.27) 



APPENDIX C: Quadratic frequency dependence 

In this appendix we extend our truncation by adding to 
eq. psp a term quadratic in the r-derivatives 

AFfe I rd^^. (C.l) 

The inverse propagator matrix involves now 



= / dxx 2 ^/l — x — 



2tt Vd 
d Vd~i 



We therefore find for the anomalous dimension 
d-K . 



(B.21) 



(B.22) 



P^A 



q^ + Vq^ + u' + 2pu" , -Sqo 



Sqo 



, q'^ + Vql+u' 



and the flow equation for V is defined by 

2A'dq'Q 



(C.2) 



dtV = vV + i-,io2dtP22\,^o. (C.3) 
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In order to evaluate eq. (jB.ip we take into account the 
modification of the propagator 



G = A-'^der^ 

( q^ + Vql+u' , Sqo 

[ -Sqo , q^ + Vq^+u' + 2pu" 



(C.4) 



with 

det = {f + Vql+u'){f + Vql + u' + 2pu") + S'^ql (C.5) 
and 

{&)n = det-2 [if + Vql + u'f - S\l] , (C.6) 
{G%2 = A-^Aet'^[{e + Vql+u' + 2pu"f-S\l], 
(G2)i2 - -{G'')2i^2A^ Act-'' Sqo{q^ + Vql+u' + pu"). 



In the definition (|C.3[1 the spacelike external momentum 
is taken at a vanishing value, q = 0. We can therefore 
replace in all propagators in eq. (jB.ip q^ fc^ and perform 
the q ' integration 



^dMq') = ^^d{k^-q'-')dt{A{k^-q''')\ 



Expanding eq. (jB.ip to second order in go we find the flow 
equation for V 



d + 2j ^ 



det 



-3 



90 



d V d + 2^ 

I - y detfc +4I/2q;,2(/,2 _^ yqf2 ^ ^ 

-5'2[fc^+u' + pM"-3y(7;)2]J. 



(C.8) 



Here det^ replaces in eq. (jC.5|) q^ ^ k^. In the disordered 
phase one finds 9tT^ = in agreement with eq. (|^ . 

In the ordered phase we evaluate eq. (jC.Sp at the mini- 
mum (u' ~ 0, u" = A) and use the integrals 



detriV^lor = T^-^/'fc'"~^"+^A„,™(«;, 5) (C.9) 



with 



and 



S 



kVV 



(C.IO) 



^n,m(w,s) = — / da;x2'"[(l + x2)(l + w + x^) + 
zvr ' 



2^ ■ svr". 



(C.ll) 



This yields 

{ - ^2,0 + 4 ((1 + ^3,1 + ^3,2) 

+ 1)^3,0-3^3,1]}. (C.12) 

The computation of 77 in appendix B remains essentially 
unchanged and one finds from eq. (|R22)) 



7^ ^ ^wXk^-^{Vk^)-^^^A2,o{w,s). 



(C.13) 



It is instructive to investigate the limiting cases s — > 
and s ^ 00. For s — the propagator (jC.4[) becomes di- 
agonal. The action has now a "relativistic" dynamic term 
involving two i9T--derivatives. This model is well understood 
and corresponds to the classical 0(2)-model in d+1 dimen- 
sions. The euclidean space-time symmetry SO{d + 1) ob- 
tains by a simple rescaling r = V^^'^t', Vd^ — d^,. The ef- 
fective action F^^o should respect this enhanced SO{d+ 1) 
symmetry. Also the flow equations would automatically re- 
spect the SO{d + 1) symmetry if we had chosen a cutoff 
function Rk consistent with this symmetry. However, our 
cutoff (IT71) violates SO{d + 1) since it only acts on d mo- 
menta q. The flow will therefore not respect SO{d + 1) 
and the full symmetry should only appear for fc ^ where 
the cutoff effects are absent. In particular, this implies 
that V{k) should approach a constant for fc ^ 0. We may 
discuss this issue in some more detail. 

Combining (|C.14p with (|C.12|) yields for s = 



dtlnV = 



1 - 



-I' 

V 



d + 2 



A2.0 



(C.14) 



1 + 1)^3,1+^3,2) } 



d + 2. 
where we define 

A = Afc''-2(V^fc2)-i/2. (C.15) 

We note that the integrals A„,m obey relations 

_d_ 

dw 
d 



-n(An 



+ l,m 



A. 



n+l,m+l 



Ar, 



-nA. 



n+l,m+l ■ 



We also employ the integral 

00 

-!- / dx{x^ + a)-^{x^+f3)-^ 



in order to compute 
1 



2(a^//3 + /3^/S) 



(C.16) 



(C.17) 



^i,o(w,0) 

^2,0(W,0) 
^3,o(?«,0) 



(i + ?« + VTT^)"\ (C.18) 

(1 + w + Vl + + W + 3%/TTw), 



16 



{1 + 110 + \/i + wy 



{5(2 + + w + 11 + llw + w^}. 
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Using eq. (|C.16[) one then obtains 



A2AW, 0) = -(1 + w + ^/l + iu)-'\2 + w), 



A^^{w, 0) + A3,2(t", 0) = -(1 + w + VTT^)-4 

o 

(5 + 2«; + i(l + + ^(1 + w;)"'/'), 



16 

{w^ - 3w; - 7 + 5(1 + _ ^ w)!/^}^ (c.l9) 



and 



1 + -j ^34 + ^3^2] =-(l + w + + 



[w^ + + 23w + 28 + 5(1 + wf^ 
-4(1 + w;)3/2 + 27(1 + u.)i/2}. 



(C.20) 



Let us consider large w where the terms ~ rj/{d-\- 2) can 
be neglected. We note that the leading term in the com- 
bination (jC.20[) ~ l/(4'u;2) cancels precisely the same term 
in ^2,0- While -q decreases for large w as Aw"^, the leading 
term in 9f In 1^ ~ \w~^ is suppressed by an additional fac- 
tor w^^ . This feature is consistent with the requirement 
V{k — 0) ^ Vq. Establishing for our cutoff the asymptotic 
constancy for V{k — > 0) for arbitrary initial conditions with 
S — ^ Q has not yet been done. We simply recall that 
any valid truncation must obey this property due to the 
50(^-1- 1) symmetry. 

For large w we find for the anomalous dimension 



d-3 



For d = 1 this yields 



dVVw 



1 



4ttVVpo 



(C.21) 



(C.22) 



and we recover the well known formula for the classical two 
dimensional 0{N) models 0,111, with K = \/Vpo. This 
is closely linked to the flow of the coupling = 1/{2k) 
in the non-abelian nonlinear cr-models [9|, |lj] , which obeys 
(TV = 2M) 



2 ^"2 4 



(C.23) 



A perturbative expansion for small or small k^^ becomes 
possible. 

In the opposite limit s —+ 00 we can neglect in eq. (jC.8[) 
all terms involving V such that (with eq. (jB.23|) and A = 



dtV 



16vd 



d+2 



d + 2 



u;Afc«^3|^l + |) / de 



(C.24) 



The r.h.s. is negative such that V is driven to positive 
values if we start with a microscopic value V{A) ~ 0. In 
this regime we find 

/VP' 



dt 



\ 52 



av 



3vd 
d 



1 - 



V 



d + 2 



u;(l + -) (1 



(C.25) 



For (approximately) constant w and A one observes two 
qualitatively different behaviors. For 775 > —1 the evolu- 
tion of 1^ ^ tends towards an infrared stable partial 
fixed point 

-2 <^v 



2(1 + 77s) ■ 



(C.26) 



In contrast, for 775 < —1 the combination s^^ increases 
fast to large values. The linear dynamic term ~ S becomes 
subdominant as s approaches zero according to 



1 3 

dts = -(1 + r]s)s + -ays 



(C.27) 



Values 775 < —1 therefore suggest a crossover from an ini- 
tial evolution where the term linear in dr characterizes the 
effective action to a "relativistic regime" where the term 
quadratic in dr dominates. If this happens and the flow 
is not stopped due to w — > 00 one expects the long dis- 
tance behavior to be governed by the relativistic model. 
For the infrared physics of the Goldstone modes the rel- 
ativistic regime applies for 775 < 0. Our findings suggest 
that this is realized for d < 3. 

Of course, once s is small eq. (jC.27|) is no longer quan- 
titatively correct since it was obtained in the limit s —^ 00. 
Also the computation of 77s has now to be performed in the 
relativistic regime. For small s we can use (jC.12[) 



dts = -(1 + 77s 



ldt\nV)s 



(C.28) 



and observe that the s-independent term in dt In V should 
be small due to the relativistic 5*0(^-1- l)-symmetry, while 
the s-dependent term is negative, reproducing qualita- 
tively eq. (jC.27|) . One concludes that the criterion for 
s(fc — > 0) ^ remains rjs < ~1, whereas the criterion for 
a relativistic Goldstone regime s <C \/w applies for rjs < 0- 



APPENDIX D: Fixed point properties for 1/ = 

In this appendix we briefly address some properties of 
the possible fixed points for w 7^ in truncations with 
a linear r-derivative, i.e. for V ^ 0. We start with the 
simplest truncation ((72)) . In terms of the variable 



Vd Xw 



d VT+- 



(D.l) 



the condition dfW ~ reads 



3 2 
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2 1 + w 2(l + 'n;)2 
3 9 



{d + 2){l + w) V (1 + w) 



0-2 ^ 0. (D.2) 



25 



Similarly, dtX = requires either A = or 



a — 2 H 7Z 777(7 



2(4-6w--w2) 



■u;(l + w)2 (d + 2)u;(l + u;)= 



fj^ = 0. (D.3) 



No fixed point exists for A = , w ^ 0. One may use a 
linear combination of eqs. (jD.2p . (|D.3|) in order to express 
a in terms of w and then solve the remaining equation for 
w numerically. Alternatively, one may numerically solve 
the flow equations ([7^ for fc - an infrared stable fixed 
point can be found easily without the need of tuning initial 
conditions. One finds a fixed point for all d < 2, as shown 
in table 1. 

As we have discussed in sect. I Villi the properties of fixed 
point (C) are strongly affected by the inclusion of the rel- 
ativistic kinetic term V. It disappears for d > 1 and 
remains for d = 1,M = 1. Nevertheless, the truncation 
= may be relevant for the initial running before a 
sizeable V is built up by the flow. It is therefore interest- 
ing to know to what extent the properties of fixed point 
(C) are robust with respect to extensions of the truncation 
which keep V = 0. For small ri/{d + 2) the terms ^ 
in eqs. (jD.2p (jD.3[) are subleading. The fixed point value 
is strongly influenced by the relative size of the con- 
tributions linear in a in eq. (jD.2[) . While the radial mode 
contributes with a negative sign, the Goldstone mode gives 
a positive contribution. This allows us to roughly estimate 
the effects of extended truncations. Adding the coupling 
7 — u'^^\pq) will enhance the weight of the radial contribu- 
tion, thus disfavoring very high values of . On the other 
hand, a contribution ~ YpQ in the wave function renor- 
malization of the radial mode will diminish its weight. We 
have investigated in appendix E the role of the coupling 7. 
While fixed point (C) persists, its location becomes rather 
unstable with respect to the order of the truncation as d 
approaches two. Only for d = 1 the fixed point seems rel- 
atively robust. 



we obtain the scaling form of the flow 
da = [2(d-l) + 37? + 27/s]7 



(E.3) 



+- 



A 



d + 2j {1 + wy 



2A^ 

1 + w 



We note that for A 7^ the flow has no solution 7 = 0. For 
fixed point (B) for d < 2 with = , A* 7^ , rj = rjs = 
one finds an IR-attractive fixed point for 7 



7* = 2A2[2(d-l) + ^A.]-' 



A 



5 - 2d 



(E.4) 



In the SSB regime we define 7 = u'-^-'(po) and use again 
eq. (|E.2[) . resulting in the flow equation 



da = [2(d-l) + 37; + 277s]7 

^2vd V \ A3 

d 



1 - 



d + 2y + 



:R, (E.5) 



with 



A^ 



and 



15 , 15 , 
— H w"*. 



(E.6) 



A, = 12 - A2w 



A3 = 6w - 12u;^ 



33 
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-w 



-w 



3 4 

■ 



(E.7) 



APPENDIX E: 
Extended truncation with six point vertex 

The neglected third derivative of u with respect to p 
contributes to 775 (cf. eq. (|B.26[) ) and to the running of 
po (cf. eq. (jA.3p V In this section we will extend the 
truncation by including u^^\ while still neglecting m^"*) and 
higher p-derivatives. Within the approximation of pointlike 
interactions (and neglecting V), the formulae for dtpo,r] 
and rjs are then complete, while the neglected term u^^"^ is 
missing in dtX as well as dtu^'^\ 

In the symmetric phase we define 7 = u*^"^-* [p = 0) and 
infer from eq. (jA.5|) the flow equation 



5(7 = 3777 



dS 



d + 2 
2A2 



(P + m2)2 \J /j2 
In terms of the dimensionless coupling 



7 



7 . 2(rf-l) 



(E.l) 



(E.2) 



For w = we recover eq. (|E.3|) . For w ^ 00 one finds for 
the leading term ~ w^^^ 



a,7 = 27ys7-^AW/2 
4d 



15-27-^-f 3 ( J- 
A2 VA2 



A2 



(E., 



While the anomalous dimension r/ is not affected by 7 
we find for 775 a correction (|B.27[) 

This contribution is positive and increases ~ w^^^ for large 
w. Combining with eq. (|E.8p yields for the flow of 7 at 
large w 

a ~ \3 1/2 



15- 15X 
A2 



(E.IO) 
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Again, for fixed A* and this yields an /i?-stable fixed 
point for and therefore for 7. In view of the fixed point 
behavior for w ^ and w 00 it may not be surprising 
that numerical solutions of the flow for d < 2 will show an 
infrared stable fixed point (w* , A* , 7*), corresponding to 
(C). 

We finally need the corrections to the flow of w and A. 
In the symmetric regime they vanish, just as the correction 
to 775. The six-point vertex does not influence the lower 
vertices. In contrast, we find in the SSB regime 



A(ftpo) - ( 1 



7 



d + 2j X (1 +w)3/2 



Po- (E.ll) 



For the evolution of A we now have to include the effect of 
the change of the location of the minimum 



dtX = dtu"{po) + jdtpo. 



(E.12) 



This yields 



3w f W7 
A2" 



(E.13) 



and 



2au 
d 



1 - 



V 



d + 2 



1 



'6w + —w^ 



A (1 + w 
' 4A2 



)5/2 



(E.14) 



We conclude that the fixed points [A) and {B) have the 
same values and A* as computed in the simple trun- 
cation 7 = 0. For fixed point {A) one has 7* = and 
the 7-direction is /i?-stable for c? > 1. For fixed point {B) 
the value of 7* is given by eq. (|E.4[) and the 7-direction 
is 7i?-stable for d < 5/2. As before, the quantum phase 
transition corresponds to (B) for d < 2 and to (A) for 
d > 2. The location of fixed point (C), however, depends 
on the truncation. The values for the extended truncation 
are shown in table 2. 



d 


w* 




7*/A2 


V 


Vs 


1 


3.22 


35.6 


0.45 


2.1 


-3.14 


1.9 


3.69 


123.5 


0.48 


2.07 


-3.99 


1.99 


3.74 


140.6 


0.49 


2.06 


-4.07 



table II: Fixed point values for (C) in truncation. 

We observe that in this truncation the fixed point comes 
close to a value where rj — 2 , i]s — —{d + 2) for which the 
fluctuation effects are relatively weak and the running of 
Po I 1 7 is therefore slow. This may well be an artefact 
of the truncation and an investigation beyond the approxi- 
mation of pointlike couplings becomes necessary. In partic- 
ular, we note that in the present truncation the expansion 
in small powers of A near d = 2 gets modified. For dtX 



and large w the terms ^ A^w^/^ do not cancel anymore if 
7 ^ A^, implying dfX ^ A once w ~ A~^. In any case, the 
large negative value of r/s indicates that the flow around 
fixed point (C) will strongly be affected by the "relativistic 
dynamic term" ~ ¥1^*8^1^, as discussed in sect. IX. 



APPENDIX F: Scaling behavior for large w 

In this appendix we investigate the flow equations (|72[) 
for large values of w. In particular, for d > 2 the flow for 
fc — !■ necessarily ends in this region since no fixed point 
with w ^ is present even in the simplest truncation. 
Including the "relativistic dynamic term" ~ V one expects 
that for all d > 1 the flow in the ordered phase obeys 
w{k ^ 0) ^ 00. We study here the simplest truncation 
and comment only briefly the extended truncations. A 
comparison of the results for the simplest truncation in this 
appendix with the results for a relativistic dynamic term in 
sect. IVIIII demonstrates the importance of the relativistic 
dynamic term. 

In the limit w — > 00 one observes for the equation (j72p 



dtw 
dtX 



w 



'^^^^"^ ^ d^d + i) 



2v 



Xid~2 + ^Xw-^/^+ ^. 

d d^{d + 2) 



The evolution of w depends on 

c 

which obeys 



^Az«V2 
2d 



2 

3w 



c- 



4(2 



d + 2 3{d + 2)w 



(F.l) 
(F.2) 

(F.3) 
(F.4) 



For large w the last term can be neglected and we obtain 
a simple closed equation for the flow of (. 

For d > 3 the combination C vanishes in the infrared 
(fc — > 0) and we obtain the leading behavior for c? > 3 



A - fc 



d-2 



X S, Po ^ S 



(F.5) 



Here the case d 
running of C 



C(fc) = 



3 is special due to the slow logarithmic 
C(fco) 



1 



(C(fco)/2) In(fcoA)' 



(F.6) 



where actually w ^ fc^^+^ , X ^ k , po ^ C'^/S- For large 
w the leading term for 775 is given by eq. ([7T|) rjs — —C 
such that S approaches a constant in the infrared for d > 3. 
For d = 3 one obtains S ^ (^'^ such that S and A vanish 
logarithmically according to eq. (|F.6p . 



5 - A - In' 



(F.7) 
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Also 77 vanishes for w —^ 00 and the infrared behavior of the 
flow simply stops, with fixed po and A. For d > 3 the flow 
is ultraviolet dominated such that = 2A(fc = O)po{k = 
0) ~ A(A)po(A) = tr, corresponding to a critical exponent 
1/ — 1/2. We recall, however, that only appears in the 
exponential decay of Gn for r 00, while 6*22 shows a 
power like decay given for 77 = by eq. (|49p . For large 
r the correlation function {(/)* {r)(t){0)) = ^{Gn + G22) is 
dominated by the "Goldstone contribution" G22 • 

The situation is different for d < 3. Now the flow of C 
exhibits a partial infrared flxed point 

C*=2(3-d). (F.8) 

The flow of w obeys near this flxed point {d ^ 2) 

dtw = (4 - 2d)w (F.9) 

and we find the asymptotic behavior 



Po 



^-1^2(3-d)^ 



(F.IO) 



For d > 2 the asymptotic value of w increases and the 
fiow always ends in the regime where w — > cx). For d < 3 
the infrared behavior of S depends on k and wc obtain for 

C = C* J w ^ 00 



(F.ll) 



According to eq. ([53)1 we note the modified relative scaling 
of time and space in the extreme infrared, 



2(d- 2). 



(F.12) 



The anomalous dimension 77 vanishes in this limit. We 
conclude that po reaches a constant value, while A and S 
vanish asymptotically 



Eq. (|F.11|) suggests that 77s increases monotonically with 
d. We may define dr such that r;s > — 1 for d > d^. Eq. 
(|F.11[) would imply dr = 2.5. For d < dr the relativistic 
kinetic term V dominates and the asymptotic equations 
for large w discussed in this appendix loose their validity. 
As we argue in sects. IVIIIl IIXI the linear kinetic term 
~ S becomes subdominant for the infrared behavior of the 
Goldstone boson physics whenever S vanishes, i.e. for rjs < 
0. We conclude from eq. (jF.lip that this happens for all 
dimensions d < 3. 

Even for d < dr only the infrared behavior in the dense 
regime is modified whereas other features continue to be 
reasonably represented by the simplest truncation. As an 
example we may consider the correlation length for the 
radial mode. For the definition of the correlation length 
we include only fluctuations with momenta q ^ > k^/^ = 
^]^^, resulting in = 2A(fci)p(fci) with w{ki) — 1. The 
behavior for tjj 3> 1 corresponds to momenta q ^ 
and does not affect the scaling ^ a^^l"^ . 



The particular scaling properties for 2 < d < 3 mainly 
affect the extreme infrared behavior of the propagator G\\. 
For g 2 — + we may replace 2A/9o ^ 2A(g)/3o with 



m = A(fci) ( ^ 



(F.14) 



We expect a similar qualitative behavior even in presence 
of the relativistic dynamic term. The value of rjs may be 
modified, however. We argue that the infrared physics for 
d < 3 is actually described by the d + 1-dimensional 0(2) 
model. This suggests r/s — d — 3 instead of 2{d — 3). This 
yields the leading part of the static propagator Gn 



r r 1 f'l' 



(F.15) 



Summarizing our simplest truncation for the ordered 
phase, we find that for d > 2 the flow drives w always 
towards large values and the condensate po or po settles at 
a constant value. For d > 3 also the interaction strength 
A and the coefficient S reach constant values, whereas for 
2 < d < 3 both A and S vanish asymptotically. For d < 2, 
in contrast, the flow drives w towards a fixed point value 
w^,. Indeed, starting with very large w eq. (|F.9[) implies 
for d < 2 a decrease of w. This continues until corrections 



w 



begin to be important. 



For the upper critical dimension d ~ 2 the situation is 
special. The coupling A becomes now a marginal coupling. 
For small enough A the running effectively stops and we can 
take a small A as a free parameter. Expanding in powers 
of A yields 



dtW — w I —2 



A 



327r VT+w 
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dt\ = (7«- 2)2(1 +7«)-5/2^, (F.16) 

One finds an infrared stable (approximate) fixed point for 
large w 



/647rY 
\3\) 



(F.17) 



As A moves slowly towards zero, increases to infinity. 
Indeed, for this fixed point the evolution of A follows 



dt 



A 



3 / A 



167r / 4 \ 167r 



(F.18) 



implying a very slow running once A(fc) ^ IBtt 




(F.19) 
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Here fcg denotes the scale where w Wt. becomes vahd. 
The effective anomalous dimension is very small 



' 2 I 167r 



On the other hand, the flow equation 
dtS = 25* , ris = -2 

implies 



(F.20) 



(F.21) 



(F.22) 



corresponding to the limit d ^ 2 of eq. (|F.11|1 . Similar to 
(IF.13P we find the asymptotic behavior 

A^AS'-fc^. (F.23) 

The order parameter is approximately constant 

fc2 (647r)2 _ (647r)2fc2 



25 9A3 18Sn\^ 



(F.24) 



More precisely, the tiny running of po for large and d = 2 
may be directly inferred from eq. (|27p (up to corrections 



dtpo 



A/Oo 



3^(l+W*)-3/2 

9AVo 



Po 



(64)37r4 



In^^- 



The evolution of po{k) stops for fc 

1 1 



po{k) = po{ko) exp 



In In # 



(F.25) 



(F.26) 



The flow of the bare order parameter po — po /A 

dtlnpo^dt\npo + v (F.27) 
is dominated by rj. For very small k < ki, 



k, 3 \ IGtt 



(F.28) 



one has 

riK\n-\kolk). (F.29) 
This implies that p{){k) vanishes logarithmically for fc ^ 

ln(fco/fc;) 



Pa{k) = po{ki) 



Hko/k) 



(F.30) 



We conclude that for d = 2 no long range order exists in 
a strict sense in this truncation. In this respect the upper 



critical dimension d 
may take the limit d 



2 is similar to d < 2. In fact, we 
2 of eq. ([5T|) and observe that with 



lim ris — 

d^2 



(F.31) 



the behavior pq ~ k'^ becomes essentially logarithmic as rj 
approaches zero for d — > 2. 

Let us compare these findings with the true evolution in 
presence of the relativistic kinetic term ^ V. Without V 
the critical dimension below which order disappears would 
be dc = 2. This is shifted to dc = 1 in presence of V. 
The presence of V reduces the disordering power of the 
Goldstone fluctuations. For a relativistic kinetic term the 
effective power counting is different - the Goldstone regime 
corresponds now to the classical model in dimension d+ 1. 

We conclude this appendix by addressing the issue of the 
scale kp associated to the density. We expect that for k ^ 
kp the density plays no role. The regime where the flow is 
essentially independent of the density effects corresponds 
to w 1. For k kp one expects a transition to a new 
qualitative regime where the density matters. This may be 
roughly associated with the Goldstone regime for w ^ 1. 

The transition between the Goldstone regime for w ^ 1 
and the linear regime for w ^ 1 typically occurs for w near 
one and we may define the scale kc where w = 1 by 



2A(fcG)po(fcG). 



(F.32) 



The detailed relation between kp and ka may depend on 
the dimension, but we expect that they are of a similar 
magnitude. For a demonstration, we discuss this issue even 
in our simplest truncation. Within the truncation (j72[) we 
have found for d > 2 that pa{k 0) = pq approaches 
a positive constant. Since po scales proportional to the 
density we can define a dimensionless quantity 



such that 



L = 2Apo " 



(F.33) 



(F.34) 



Since L is evaluated for k — kc we expect i to be a con- 
stant that is neither extremely small nor large. This implies 
ka = VLpI^"^ ~ kp. 

For d < 2 the issue is more involved. Besides w the 
crossover to a relativistic kinetic term plays a role. Within 
the truncation (j72p one may wonder where the scale set 
by the density appears in a situation for which the cou- 
plings w and A fiow to their fixed point values irrespective 
of their microphysical values. If all couplings are irrelevant 
the information about the density would be lost in Ffe^o- 
Expressed in terms of renormalized fields the effective ac- 
tion should therefore contain at least one parameter that 
is not determined by the fixed point. Such a parameter 
is given by S and the presence of a scale can therefore be 
encoded in S, despite the fact that S is dimensionless. In- 
deed, the nontrivial scaling with 77s 7^ implies the generic 
form 



S — >S'n 



(F.35) 
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and a momentum scale appears by dimensional transmu- 
tation in the form of fcp . Since the size of S determines the 
scale where V will start to dominate, it is plausible that 
ko can be related to the final value of po and therefore to 
n for d > 1. For d < 1 the condensate contribution to the 
density po vanishes and the properties of G{q) for q 
play a decisive role. 

We have avoided these subtle points by choosing a fixed 
definition (|62p for kp- The price to pay is that the transi- 
tion between the qualitatively different behavior for the 
dense and the dilute regime occurs not necessarily for 
kph ~ kp, but may involve a nontrivial proportionality 
constant. 



APPENDIX G: Symmetries and thermodynamic 
relations 



In this appendix we summarize symmetries and Ward 
identities [1], Q on the level of an appropriate trun- 
cation of the effective action Let us consider, for 
real time (Minkowski signature), the classical action Sm 
in presence of local sources 



S 



M 



X 



X 



idt+cr + {V ~ iAf 

v{x'x)+rx + ix*}- (G.l) 



Here (j{x) and A{x) are real source fields, while j{x) is 
complex. For a homogeneous setting the physical values 
will be <j{x) — ^,A(x) = Q,j{x) = 0. The action is invari- 
ant under time- and space-translations, rotations, parity 
refiections Xk — > —Xk, Ak —f —A^^ as well as time reflec- 
tion t — > —t, p —f p* ^ j ^ j* ■ It is real and exhibits a 
local [/(l)-gauge symmetry associated to conserved parti- 
cle number 



A^{x) A^{x)+df,d{x) 



(G.2) 



where A^ — (u, A), — {dt,W). For yl = the action 
remains invariant under Galilei transformations 

X{t,x) ^ e'fx{t,x~2Pt), j{t,x) ^ e'fj{t,x-2Pt), 
f = Px-P^t, a{t,x) ^ a{t,x-2Pt). (G.3) 



In eq. (jG.ip we have scaled our units such that effectively 
2Mb = 1, with Mb the mass of the particles, such that 
the velocity obeys v = P/Mb=2P. Thus for a plane wave 
with uj = the Galilei transformation indeed produces 
the appropriate shifts in momentum and energy, = p + 
P, J = {p + Pf, i.e. 



J{p' x-uj't) 



(G.4) 



We note that Galilei transformations and x-dependent lo- 
cal U{1) transformations are not compatible. In the follow- 
ing we take A = Q and restrict the C/(l)-transformations 



to transformation parameters z?(t) depending on time but 
not on the space coordinates. 

Including the fluctuation effects yields the effective ac- 
tion r[<^, ct], with ^{x) — {x{x)) evaluated for a corre- 
sponding source j{x). In absence of anomalies from the 
functional measure, and for T = 0, the effective action has 
the same symmetries as the classical action, where p trans- 
forms in the same way as x- This extends to the average 
action if the cutoff Rk is consistent with the symme- 
tries. For an investigation of the long distance behavior 
we expand the effective action in the derivatives dt and A. 
Consistency with the symmetries requires 

r = j^^[Uo{p)^\z[p)[p*[idt + a + A]p + c.c] 

-'\y{p)[^* [idt + <t + aYp + c.c} 



p-2idtV{(p*Vp~pVp*) 



+2V{Vp*Ap + VpAp*) - A{p*Ap + pAp*)^ 

-\y{p)pAp + Ac]. (G.5) 

Here we use p = p*p and AC contains further gradient 
terms like [VpVpY or V pV{(p*{idt + a -\- A)(p). In par- 
ticular, Galilei symmetry implies that all time derivatives 
appear only in the combination = dt — ic — iA, and ad- 
ditional gradients act only on the invariants p and p*!)^^. 

We next perform an analytic continuation to the eu- 
clidean effective action by replacing —idt dr in eq. 
(jG.Sp . (We have introduced an overall minus sign in the 
transition from Si\i to F, such that the standard conven- 
tions for the classical action in Minkowski spacetime match 
with the standard conventions for the euclidean effective 
action.) This allows for an extension to thermal equilib- 
rium with T 7^ 0. For T ^ additional terms can appear 
in the effective action. First, Galilei invariance is broken 
since the heat bath singles out a particular reference frame. 
This permits the appearance of additional terms involving 
gradients, as p* Aip. Second, a continuation of the local 
[/(I) invariance with real functions i?(t) requires complex 
a, transforming as cr — > a + idr'd. However, the real part of 
a is invariant under local transformations. In consequence, 
an arbitrary dependence of U, Z, V on a+a* becomes possi- 
ble. In practice, only the global C/(l) symmetry is relevant 
for T =/= 0, since the local transformations only constrain 
possible couplings of the imaginary part of a which play 
no physical role. The Ward identities related to the local 
C/(l) symmetry and the Galilei transformations therefore 
only restrict the limiting behavior of F for T ^ 0. We 
emphasize that all Ward identities are automatically im- 
plemented by the invariant form of F in eq. (jG.Sp . 

The effective potential in presence of a nonzero chemical 
potential a reads 



Uip, a) = Uoip) ~ Z{p)pa - V(p)pa^ 



(G.6) 



The order parameter is given by the condensate density 
= POi which corresponds to the minimum of J7 at a 
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fixed value of <j 
dU 

— = , apC/o,,„ = adpiZp + aVp)i,,. (G.7) 

The total particle density is related to the cr-derivative of 
U at po 

dU 

n - —g^ipo) - [Z{po) + 2aV{po)]po. (G.8) 
In terms of renormalized fields 

cj) = A^/^, p = Ap, Z{pq) + 2aVif5o) (G.9) 
one finds 

n^po. (G.IO) 

Thus a nonvanishing density requires for T = a nonvan- 
ishing renormalized order parameter. (We recall that there 
are p- independent contributions to U for T ^ 0. For exam- 
ple, a piece AU = —riTcr contributes to n a piece ut that 
does not vanish for pQ = 0.) This observation is partic- 
ularly interesting for d = 1 where the consendate density 
Po vanishes in the infinite volume limit. Indeed, we find a 
diverging A as the IR-cutoff k runs to zero. 

We may also compute thermodynamic susceptibilities 
like the response of n to a change in the chemical potential 

OCT |T=o OCT oa 

S^A\l + ^^^^^ uw^ \ g-^^ 

y a In po ) 

Differentiating the minimum condition (|G.7p yields 



dpQ _ 1 dn{po,c7) - _ d'^U{p,a) 
da A dpQ 1^ ' 9p2 

where 

dn d{Z{po)po) , „ d{V{po)po) 

— zcr- 



(G.12) 



9po \a 9po 
and therefore 



dpo 



S, (G.13) 



dpo^S dn _ 2Vipo)n 

da A ' da\T=o A ^ A ■ ^ ' 



The pressure 



P = -U{po) 



(G.15) 



is normalized such that it vanishes for T = n = 0, i.e. 
Uo{p = 0) = 0. Its response to a change of the chemical 
potential obeys 



dp dU dU dpo 

da\T da\po dpo\^^p„ da 



(G.16) 



In our units {2Mb = 1) the energy density obeys e = n/2 
and one obtains for the macroscopic sound velocity 



dp 2 ^ 



* de \T da\T \da\T 

For T = one finds 

-2 ^ £(Po) _^dpo ^ Vipol ^ 
'^'^ A 2po da 



A 



2poX 



(G.17) 



(G.18) 



In presence of spontaneous symmetry breaking we are 
interested in the propagator for small fluctuations around 
the expectation value (po that we take here to be real. We 
decompose (p{x) = <po + '^{Sfiix) + i'p2{x)) . The in- 
verse propagator matrix is encoded in the part of F that is 
quadratic in Sipi and (p2- Retaining terms with up to two 
derivatives one may parameterize 



(G.19) 



F2 = J ^^{Sipidr(p2 - (P2dTS(pi) 

-^{(P2d'^(p2 + BrSipid'^Sifi) 

1 JTl^ 

--{1P2A1P2 + CrScpiASipi) + -^Sipj 



Here we have defined the renormalized field ip such that 
the term ~ ip2Aip2 has coefficient one. For T = we 
may extract the different couplings by expanding eq. (jG.5|) 
(with A = 



F2 = 



Z{po) + 2aV{po) + poidpZ + 2adpV)ip„ 



{6ipidrip2 - (p2drSipi) 

-\v{po)(p2dl(p2 



[V{pq) + po(2 + dpV + Yt+ pdpYt)\p^]5<pidl5^i 



[Z{pq) + 2aV{pQ\Lp2ALp2 

Z{po) + 2aV{po) 
+po{2dpZ + ^adpV + Y + pdpY)\p„] {.S^iA5<pi) 
+\p^&V>l]- (G.20) 



The renormalized field is the same as in eq. (|G.9p and we 
identify 



5* 
V 

Br 



^ ^ dln{Z{po) + 2aV{po)) ^ S_ 
9 In Po A' 

Vjpo) 
A ' 

^ ^ g ^lnV'(po) ^ Yt{Po)Pa ( ^ 
d\npo 

25-l + ^%l^(l 
A 

2Xpo/A = 2Xpo , A 



dlnYtjpo 
iViPo) V 91npo 

ainr(p)^ 



a In Po 



2^ 



(G.21) 
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We will find that S vanishes in the infinite volume limit 
for d < 3. The solution of eq. (|G.21[) implies that A{p) 
diverges for p ^ 



A{p) = Zip) 



2aV{p) - ^. 

P 



(G.22) 



For d > 1 one finds a nonzero po such that A — A{pq) 
remains finite. For d = 1, however, po ~^ and A diverges 
in the infinite volume limit. Since V remains finite also V 
must diverge ~ p"^ in this case. 

Comparing the definitions (jG.lip and (|G.2ip one has 



c-' = V 



2Xpo 



(G.23) 



For S* — > this implies for the macroscopic sound velocity 

(G.24) 



In summary, Galilei and gauge symmetry relate for 
T = the properties of the inverse propagator, like the mi- 
croscopic sound velocity v, the superfiuid density ns = Apo 
(as defined by the stiffness with respect to phase changes), 



or the term linear in the frequency S, to macroscopic 
thermodynamic quantities. Using 



dn/da = 2Vn + S^/X = 2n/cl 



(G.25) 



we can replace the cr-derivatives by n-derivatives or deriva- 
tives with respect to the volume fi^ {dlnn/dlnfld = — 1 
for fixed particle number N) 



din 'I 



2 dlnrir 2 91nnc 



da 9 In 7 



2^ d In ild 



(G.26) 



This gives a direct physical interpretation of the running 
renormalized couplings evaluated for fc = 



po = n ^ Us , po ^ Uc , A 1 = — = 



S 2naini 



n 



A 91nn ' 



9 Inn 



(G.27) 
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